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ABSTRACT 


The S-Matrix is defined as an the limit of the Time evolution operator U(te,ti) between 
asymptotic initial and final Hilbert /Fock spaces for a time-dependent scattering process, i.e. 
between too and t—-oo. There unitarity encodes conservation of probabilities over time. 
The goal of this thesis is to do a thorough and complete analysis of this useful quantity in 
physics namely the S-Matrix, learn its construction specifically in the path integral formu- 
lation, and apply it in the forefront areas of Quantum gravity and String theory.The most 
trivial construction of S-Matrix can be done as an infinite perturbative series(Dyson Series) 
in the Canonical procedure of QFT.An analysis up to first order is done by introducing in 
-states and out-states and a Born approximated result is obtained. Further both the inter- 
nal and external symmetries of the S-matrix are studied and their implications are drawn. 
An analytical proof of the unitarity of S-matrix is presented and further its implications are 
studied. A brief introduction to the Path integral formulations of QFT is done in order to 
construct the S-Matrix in the Path integral formulation, for this purpose various generat- 
ing functionals are analysed. This is because greens functions are nothing but functional 
derivatives of the generating functionals also called partition functionals.Connected and am- 
putated greens functionals and their generators are studied. As expected the S matrix in the 
Path integral formulation turned out to be a functional rather than a series or operator as 
in canonical QFT. After studying a general prescription to construct the S-matrix in Path 
integral formulation its unitarity is studied.This is very important if we want to construct 
any theory of quantum gravity then the S-matrix constructed for asymptotic interactions of 
the theory must be Unitary in order to preserve probability and respect causal conditions. 
At last introductory aspects of Quantum gravity and various approaches to quantize gravity 
is studied. The S-Matrix of Quantum gravity is also studied by Canonical Quantization of the 
Gauge Field.Introductory aspects String theory is studied via the Path integral formulation 
(Polyakov Integral) as it seems to be the natural language to construct a theory of Quantum 
gravity.Finally an attempt to construct S-matrix of these theories are done where the last 
sections are tributed to analyse the Unitarity of the S- matrix of the current string theory 
is done.As we know that the S-Matrix of these theories also needs to be consistent with the 
Unitarity criterion. However the Non Unitarity also seems to creep in some areas of standard 
model as well as 2D string theory,but that is beyond the scope of this work.Our sole aim is to 
provide a detailed analysis of the S-Matrix,the physics that comes with it,its uses, Unitarity 


and its implications. 
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1.1 Introduction 


The main goal of this thesis is to do a complete analysis of one of the most important 
quantities in high energy physics i.e the S-Matrix.The analysis will be comprehensive 
thorough and self sufficient.We start with a formal introduction of the S-Matrix and 
discuss the major differences in two major formulations of modern Quantum field the- 
ory(QFT)namely, the Canonical formalism vs the Path integral formulation.Although 
the former i,e the Canonical formalism is used often by theorists as it seems most nat- 
ural way to quantize a theory by invoking canonical commutation relation and seems 
sufficient. However we have learnt that although this approach is very powerful it seems 
to fail beyond the standard model, and is miserable enough when applied to construct 
theories of Quantum gravity.Path integral formalism on the other hand seems to be 
more powerful in approaching the modern problems in String theory and Quantum 
gravity.This is expected because the path integral formalism has its roots in one of 
the most fundamental laws of nature i,e the extremization of the action or the action 
principle itself. This is what has been noticed and applied by many giants,Feynman 
and Hawking to name a few.Path integral formulation can be easily generalised to the 
theory of Quantum gravity and String theory ,although the mathematics seems to be 


very tedious but the results seems promising. 


1.1. Introduction 





Our study starts with the analysis of the S-matrix in the canonical QFT in chapter 2 
where we study its properties ,symmetries,its uses and its general structure. We will 
try to visualise its action in the matrix form as well as in physical space, and we finally 
obtain its expression in Perturbative form and an Analytical expression of Ist order.At 
the end of the chapter we discuss one of most sought after properties that the S-Matrix 
of any theory must have i,e its Unitarity and study its implications in the canonical 
QFT and Scattering theory. 


We introduce the machinery of Path integrals in chapter 3 where in first few sec- 
tions we discuss the general aspects and workings of the Path integrals in the regime 
of NRQM(Non Reltivistic Quantum Mechanics ,Just to get and idea and make the 
basics clear and then we move to the Path integral formulation of field theories.In 
the field theory Path integrals we will study important tool namely the Generating 
functionals which will help us in the construction of the S matrix in the Path integral 
formalism. We will also study the generalised n-point correlators and path integral 


formulation of n-point Green’s Function. 


Having developed all the tools ,we now start constructing the S-Matrix in the path 
integral formalism of QFT in chapter 4. We will extensively used the tools developed in 
the previous chapters. We will also dedicate one subsection for analysing the Unitarity 
of the S-matrix in the Path integral formalism. This is important because unlike the 
Canonical formalism the S matrix here is a complex functional instead of an operator 
that we are used to in the canonical formalism. It is important to keep in mind that 
unlike in Canonical formalism where most of the objects we deal with are some kind of 
operators we don’t have any concept of operators in the path integral formalism.Most 


of the dynamical variables are functionals in the Path integral formalism 


Finally in the last chapter we will apply all that we have learnt in the previous chapters 
and study the Path integral formalism of the Quantum gravirty and the String theory. 
As already mentioned this will be done in the Path integral formalism because it seems 
to be the natural language. At last we will try to construct S-matrix in both Quantum 
gravity and String theory by drawing analogies from chapter 4. It is worth mentioning 
that the construction will be somewhat unrigorous and is just used to demonstrate the 
importance of the S- matrix. We will conclude the thesis by asking ourselves why is 
the Unitarity principle so important in physics and in particular what the Unitarity of 


the S matrix tells us about the theory, as we know the mordern String theory came 


1.2. Canonical Vs Path Integral formalism 





into existence by studying the S matrix theory in. 70’s so something deeper must be 


going on. 


1.2 Canonical Vs Path Integral formalism 


Historically There have been two different approaches to any Quantum Field Theory 
in the past. The first was developed primarily by Paul M Dirac in order to deduce 
Quantum Electrodynamics in 1927. Canonical Quantization is the term given to this 
method ”the procedure By advocating fields as operators and introducing quantization 
via canonical uncertainty relations, this formulation takes a more direct approach to 
QFT. This is done in order to maintain maximal symmetries in the theory. As one can 
see, the term ”canonical”’ denotes that we operate along the lines of the Hamiltonian 
formulation, which means that we already have a Hamiltonian (H) and the system dy- 
namics are generated using canonical Poisson brackets/commutators. Although more 
general commutation connections are introduced in quantizing the theory of fields, this 


aids in the preservation of not only exterior but also internal symmetries. 


The second method is referred to as the “Path integral method ” in which Rather 
than quantizing the canonical variables q and p, and hence the canonical relations, we 
quantify a more fundamental quantity of every dynamical system, namely the action [S] 
itself.We treat particles (QM) or fields (QFT) as if they were simultaneously travelling 
or evolving through all possible paths in this formulation. Despite the fact that the 
mathematics is more complicated, this formulation is more fundamental because it is 


based on the fundamental principles of dynamics itself. 


In the Canonical Quantization formulation a quantity called the “Scattering Matrix” (S- 
Matrix) plays a very fundamental role. Although it does not introduce any new physics 
to the formulation, it does provide additional insights into scattering phenomena and 
compliments the computations done via the Feynman diagrams. The “Unitarity” of 
this S-Matrix provides the mathematical basis for the probability conservation in all 


of Quantum field theories. 


However in the Path integral Formulation, on the other hand, one does not explicitly 
perform a comprehensive analysis of the S-Matrix.It’s a bit of a grey area, especially 


when it comes to its Unitarity condition.However, in order to preserve probability 
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(which is a necessary criteria for any formulation/ any physical theory), the sum of 
probability over all paths/field configurations must equal unity.This can be accom- 
plished by creating a correct S matrix and demonstrating the Unitarity condition for 


it. This will be the main area of focus of this thesis. 


1.3 Scattering Matrix(S-Matrix) 


The S-matrix is defined as an the limit of the Time evolution operator U(ts,t;) between 
asymptotic initial and final Hilbert spaces for a time-dependent scattering process, i.e. 
between too and t—-oo.There unitarity encodes conservation of probabilities over 
time.More formally, in the context of QFT, the S-matrix is defined as the unitary 
matrix connecting sets of asymptotically free particle states (the in-states and the out- 
states) in the Hilbert space of physical states. A multi-particle state is said to be free 
(non-interacting) if it transforms under Lorentz transformations as a tensor product, 
or direct product in physics. Asymptotically free then means that the state has this 


appearance in either the distant past or the distant future. 


1.3.1 The Power of the S -Matrix 
In particle accelerators, physicists measure the following quantities: 


e cross sections of scattering processes 

e masses of stable particles 

e magnetic moments of elementary particles 

e life of unstable particles 
The task is to theoretically predict the measured values. To this end, physicists mainly 
use two concepts, namely, 

e the S -matrix and 


e n-point correlation functions/Greens Functions 


4 
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Roughly speaking, the two concepts are equivalent. The S -matrix encodes the tran- 
sition probabilities of scattering processes. In fact, the S -matrix knows all about 


elementary particles. 


e The S -matrix enables us to compute the cross sections of scattering processes. 
e The singularities of the S' -matrix refer to the masses of stable particles. 


e Analytic continuation of the S -matrix in momentum space tells us the masses 


of unstable particles (also called resonances) and their lifetimes. 


e The S -matrix allows us to compute the full correlation functions by the magic 
Gell-Mann-Low formula. This way, for example, we get so-called form factors 
which yield the magnetic moments of elementary particles (e.g., the anomalous 


magnetic moment of the electron and the muon). 


While the S-matrix may be constructed for any quantum mechanical scattering pro- 
cess, we’re mainly interested in the asymptotically free particles scattering and in- 
teraction scenario. We'll address this using the field theory formulation, although a 
brief Introduction will be given in terms of the ordinary Quantum mechanics so as to 
draw analogies and provide the reader with a better insights and understanding of the 
underlying steps.We begin with a brief overview of the time evolution operator in con- 
ventional non relativistic quantum mechanics, then move on to the interaction picture, 
and ultimately to the time evolution operator for fields using the path integral method. 
Finally, we’ll set the needed restrictions on a suitable S-Matrix and then establish its 
unitarity for any general situation. The goal of this study is to provide yet another 
evidence that the S matrix is unitary even when using the Path integral formulation. 
These are the reasons for this type of action. For starters, the unitarity condition is 
one of the most basic prerequisites for a theory to be considered physically acceptable. 
Second, it’s preferable to have as unambiguous a proof scheme as feasible, with the 
process ensuring unitarity easily identified. Third, given the wide range of Lagrangians 
seen in modern quantum field theories one would like to know what criteria, to put on 


the regularisation procedure, etc., that may guarantee that the S matrix is unitary. 


CHAPTER 2 


S-Matrix analysis in Canonical QFT 
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2.1. Time Evolution 





2.1 Time Evolution 


In general, physical systems are dynamic, meaning they change over time. Newton’s 
laws of motion explain the dynamics of classical mechanical systems, while Maxwell’s 
equations determine the dynamics of the classical electromagnetic field. Although 
some unification is provided in terms of Hamilton’s principle, a universal manner of 
articulating the laws of classical dynamics, either Newton’s equations or Maxwell’s 
equations constitute laws of physics whose form is specific to the kind of system being 
discussed. We need a way to explain the dynamical properties of quantum systems, 
i.e. given the state U(0) of a quantum system at some initial time t = 0, we need a 
physical rule that tells us what the state will be at some later time t, i.e. a quantum 
law of evolution. Given that a state vector is a repository for all known information 
about a system, we need a general physical law to explain how this information evolves 
over time in response to the physical circumstances in which the system of interest 
finds itself. It’s to be expected that the details of this law will differ from system 
to system, but it turns out that the law of evolution may be described in a fashion 
that applies to all physical systems, similar to how Hamilton’s principle allows for the 
concise expression of classical dynamical laws. Schrodinger’s equation is, of course, the 
quantum law of evolution.Time evolution refers to the change in state caused by the 
passage of time in systems having internal states (also called stateful systems). Time 
does not have to be a continuous quantity in this formulation; it may be discrete or even 
finite. The laws of classical mechanics control the temporal development of a group of 
rigid bodies in classical physics. These principles, in their most basic form, describe 
the connection between forces acting on bodies and their acceleration as determined 
by Newton’s equations of motion. Hamiltonian or Lagrangian mechanics may also be 


used to explain these concepts in a more abstract way. 


Other stateful systems may benefit from the notion of temporal evolution. Stateful sys- 
tems are often described in terms of both states and observable values. Time evolution 
may also refer to the change in observable quantities in such systems. This is especially 
important in quantum physics, where the Schrodinger and Heisenberg pictures of time 
evolution are (largely) identical.Time, unlike position, is not an observable quantity in 
quantum physics. There is no Hermitian operator with eigenvalues equal to the sys- 


tem’s time, and time is just a parameter, not a measurable entity. So, according to the 
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relativity theory’s teachings, time and position are not on a same footing. The equality 
is restored in relativistic quantum field theories by degrading the position down to the 
parameter level.In short we first have to understand the process of time evolution and 
for that it is necessary to introduce the mathematical formalism of the “Time evolution 


A 


operator” [U(t,,ti)] 


2.2 Time Evolution operator 


In Quantum Mechanics, we know that every system has many quantum states. As a 
passage of time new states are evolved. What we wish to do is to find a relation between 
these initial states and the evolved states. For this purpose we use a linear operator 
which evolves the states over time which is famously called the Time Evolution 
Operator. Let us see the various representations of this time evolution operator in 
different pictures and understand how it plays a fundamental role in ensuring the 
Unitarity of our final S Matrix.We start with the Schrodinger picture with a free non 
interacting Hamiltonian and subsequently move towards a more general representation 


in the Interaction Picture 


2.2.1 Schrodinger Picture with free Non Interacting Hamilto- 


nian 


Let’s start at the beginning and get the equation of motion for a system in general 
quantum state W(t) and its time evolution through the time propagator. We’re looking 
for quantum motion equations that are analogous to Newton’s — or, more precisely, 
Hamilton’s — classical motion equations. How does this state vary with time if we 
know it at time V(t;) How do we calculate V(t;)for a future time t; > t;?In this case, 
we'll rely on our intuition, which is primarily based on classical mechanics. We will 
not explicitly demonstrate the spatial dependency of the state in the following section 


to keep notation to a minimum. 


A major postulate of Quantum mechanics also known as the famous Schrédinger equa- 


tion states that for any system in a general quantum mechanical state evolves according 


2.2. Time Evolution operator 





to the equation 
dW (t) 
dt 


The most general solution to the above equation is given as 


ih = HWV(t) (2.1) 


|W (t,)) = eM 4/ DW E,)) (2.2) 


where the exponential term 
e tH (tet) = U (ts, ti) 


is known as the quantum mechanical time evolution operator.In general for a quan- 


tum particle travelling from a point (xq, tq) to another point (2, t,) the matrix element 


ta) = (a te) Q3) 


It is worth mentioning that We start by assuming causality: |W(t;) > precedes and 


U of the time evolution operator U is given by 


U (2p, te; Dasha) en i(te—ta)H/h 


Ul tyte tats) = (1 














determines |W(t;) > , which is crucial for deriving a deterministic equation of motion. 


Also, as usual, we assume time is a continuous variable: 


In case of Quantum field theory the matrix element of time evolution operator ana- 
logically known as the Transition Operator for a field (consider a scalar field ®) 


evolving from an initial field configuration®; to a final Configuration ®, is given by 


(a, 


where H = J @xH is the total Hamiltonian of the system and H is the Hamiltonian 
density. 


evilts—ti) A 








®,) (2.4) 


2.2. Time Evolution operator 





2.2.2 Properties of time evolution operator /Transition opera- 


tor(For time independent non interacting hamiltonians) 


1. Unitarity: This is the most important criterion for any time evolution operator. 
It says that 
UUt=1 (2.5) 


This criterion is must for giving correct probabilistic interpretation to any wave functionV 
or any field ®.The conservation of probability in quantum mechanics and QFT comes 
from this Unitarity condition.As we know that sum of all probabilities should be equal 


to unity. 


2. Time ordering/Causality: Another important criterion for any time evolu- 
tion operator is the ”time order” that is the order in which it acts on the concerned 
state/ field.This is important in order to prevent any theory from breaking the ” causal- 
ity” .Hence for any state/field ® evolving from timet, to time t,, where the chronological 


order is from time 


The evolution operator U(t,,t,) between the states ®(t,) and ®(t,) must follow the 


semi group property given by 
U (tn, t1) = U (tn, th).U (th, t1) (2.6) 


Note that this property directy comes from the Hermitian nature of the Hamiltonian 
operator H. 

H = Ht 
This is guaranteed by the fact that the Hamiltonian is Independent of time.(Atleast in 


the above section. 
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2.2.3. Interaction Picture and Time evolution 


Now we shift our attention to a more general case of time dependent hamiltonians, 
which indded is the case for many Interaction theories. Suppose that we have a Hamil- 


tonian H, which we split into a free part H, and an interacting part H7, i.e. 
H=H,+4+ 4H, 
We will work in the interaction picture here and assume units such that the reduced 


Planck constant h is 1. 


In the interaction picture, the evolution operator U defined by the equation 


W(t) = UE, to) V(to) (2.7) 


is called the Dyson operator. We have 
UE) =f 


U(t, to) = U(t,t1)U (tr, to) 
U-*(t, to) = U(to, t) 
Now we differentiate the equation (2.7) both sides with respect to t and put the initial 


condition i.e 


dU(t) — dU(t, to) 











= L U(t, 2. 
di di V(to) + UE, to) di (2.8) 
dWv(t)  dU(t,to 
—x 2. 
eC) (2.9) 
In interaction picture 
Uv 
no = H,(t) V(t), hol (2.10) 
a = H,(t)V(t) = —iH;(t)V(t) (2.11) 
then the equation becomes 
t 
wt fo) _ re) W(t) (2.12) 


11 


2.2. Time Evolution operator 





aU (t, to) 
dt 


Equation 2.13 is commonly known as Tomonaga—Schwinger equation integrating 


= H,(t)U(E, to) V (to) (2.13) 


the above equation from t; to ts 
ty 
i) dU (t, to) = [inoue to) dt 
ti 
CR nee ee i H,y(t,)U (t,t) dty (2.14) 
put t; = to and ts =t , since U(to, to) = 1 
Ulipto) = UG) = 1 =i f Hi(t)U (tt) dt, (2.15) 


which leads to Neumann series and the equation becomes 


t 
U(t, to) _ dt,Hy(t,) + car fan f dtyH;(t,)Hy(tg)... 
(2.16) 
(—7) oe an, f dto.. 7 dt nHy( (ty ) Hy (te). Ay (t n) +... 


since t; > tg > ...tn, fields are time-ordered and it is convenient to introduce an operator 


T called ” Time-ordering” operator which is given as 


U(t,to) = (-#) fan fe dty.. J “Gar (41) Hr(t2)..-Hr(tn)] 


By making the integration simpler we can write it as 


t ty tn—-1 
Sx — / dt, / dto... ‘ dt, K (ty, ta, saith) 
to to to 


Let us assume that K is symmetric in its arguments and define 


t t t 
Tn — / dt, / dty... / dt, K (ty, ta, sabes) 
to to to 


so we can say the region of integration can be broken into n! sub-regions by t, > tg > 


. > th, te >t) >... > tn. The integral in each of these sub-regions is the same and 
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equal to S,, because of the symmetry of K, hence we can write , 


By using the above identity ,we can write the previous integral as 


n= 





= ag\n t t t 
— [uf dts. f dtyT Hy(t1) Hy (t2)...H1(tn) (2.17) 
: to to to 


summing up all the terms,therefore we get the Dyson’s theorem for the Dyson series 


as 
U (t,o) = S> Un (ts to) = Te MoM (2.18) 
n=0 
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2.3 S -Matrix and Scattering Theory in Canonical 
QFT 


2.3.1 S-Matrix as the limit of Time evolution operator (Perturbative 


Expansion form) 


The direct way to define the S-matrix start with consideration of Interaction picture. 
Let us consider the Hamiltonian H to be split into the part Ho and the interaction 
H,;,H = H) + H;. As we know the operators behave as free field operators in this 


picture and the state vectors have dynamics according to the interaction Hy. 


Let |(t)) be a state that has evolved from a free initial state |®;). Then the S-matrix 


element is defined as the projection of this state on the final state |®,). Therefore, 
Spi = lim (®/|W(t)) = (/1S1,) (2.19) 


where S' is the S-operator. The great advantage of this definition as that the time- 


evolution operator U evolving a state in the interaction picture. 
U(t, to) = Te fio art) 
where 7 denotes the time-ordered product. Expressed in this operator, 


tg7+00 t1 4-00 


from which we can write S = U(co, —oo) expanding using the knowledge about U gives 


a Dyson series , 





= : “f dee. ae dt T[Hy(ts)... Hy (ty)] 


n= 


or ,in terms of Hamiltonian Density, 


ox ye — is diya: ie dx\T [H1(t1)..-Hi(tn)| 
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Being a special type of time-evolution operator, S is unitary. For any initial state and 


any final state one finds 





det... i : teh (Halt) Halt) ®s) (2.21) 





= (0,|5|0,) = (8 


2.3.2 The S Operator 


From Eq 2.21 , we can find any given S; component of the S Matrix, such as the f = 2 








final state, 
(f = 2||F) = (1 f=2 de Sal fay = 2| (Sul f = 1) + Soil f = 2) + S3ilf = 3) + ..) 
Sl, 
(Soper li +Soper .2i t+...) |t)=Soper |) 
= (f = 2|So;| f = 2) + ( many zero terms ) = So; = (f = 2|Soper | 4) - 


Hence In general, 
Spi = (Ff |Soper |) 


This operator Sandwiched between the initial and final states is called the S-Operator 
2.3.3 Expanding Soper 
Soper Can be expanded (in a Taylor series like e* = 1+ 2+ 27/2!4 2°/3!+...) as 
if. ! HE (t)at I 
ern (tet) = e- tare =![-— if A; (t1) dt; 


-ar((Lmn) (f meoe)} 


(—7) 
= = - = i T { Hj (t1) Hj (to)... H} (tn) } dtydta... dt 
n=0 : bi ti 


(2.22) 





If the H/ above were numeric functions of time, it wouldn’t really matter what order, 
with respect to the t,, we carry out the integrations in. However, since they are 


comprised of operators that act on a ket state to their right, , we have to be sure 
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that at each point in the integration, the time wise earliest operators are acting first. 
The operation of the operators must occur in the order they occur in nature, in time 
sequence. So, the earliest operator must operate on the ket first, i.e., it must be furthest 
to the right. The next earliest to operate must be second from the right, etc. This 
means that at each point in the n dimensional space where each axis is a different t,,, 
the H/ dependent on the earliest time of the ¢,, should act first, the H} dependent on 
the next earliest of the t,, should act next, etc. See Figure below for a two dimensional 


illustration of this, where we are integrating over the t, — ft» plane. 










t.<t, 





Integrate in this order t< t, 


Hy (t,)H; (t,) at, at, 


BN 


over t,<t, region 





Integrate in this order 


FAA Bit), (t,at dt, 


over f, < t, region 


f, 


Figure 2.1: Time-ordering of integrand in t, space 


We indicate this ordering with the time ordering operator symbol 7 in the integrands 
. The order of the integrand operators is rearranged as we integrate over all t,, dimen- 


sions, such that the operators are time ordered at every point (t), ta,...tn) 


Taking our integration limits to infinity in both space and time, equation 2.23 is turned 
into the time ordered infinite spacetime Soper , i-e., the Dyson expansion of the S 
operator, 


S = Soper (ty — 00,t; 4 —00) = ete Hide 
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Which on expansion gives 





pal [ or i‘: T {HE (ty) HE (ta)... HE (ta)} dtedtan... dtr, 


n! 


SPECIAL NOTE: Above equation is the usual definition for the term ”S operator” 
as used in QFT. Infinite time and space, along with time ordering, are assumed for 
the more general case of arbitrary time and space limits, in what we have been calling 
Soper - 90 a big difference between Soper and S is that time integration and volume of 
interaction for the latter are always infinite. From now on, when we use the term ”S 


operator”, unless otherwise stated. 


2-4 NOTE: Similarly, the S matrix elements S$; in the literature usually correspond to 
the values for infinite time duration and infinite volume of the interaction. For us, they 
could represent either case, of finite time/volume or infinite time/volume, but from here 


on, unless otherwise specified, they will correspond to the infinite time/volume case. 


ie 1 fore) fore) ; a - 
S= i “if Meaty ff TH co (ed }atrdina += 
so ———__J_—_—_— nS Se 
gQ) S(2) 


2.3.4 Ways to Visualize the S Operator 


2.3.4.1 Matrix form 


In a quantum scattering process different interactions have different probabilities to 
occur. So to keep track of each and every possible interaction, or to keep track of the 
individual probabilities for each to occur we use the machinery of S-Matrix shorthand 
for Scattering Matrix (since each of these interactions is considered a scattering 


process). 


Consider a column vector for which each component represents a different initial (in- 


coming) quantum eigenstate (typically multiparticle), as on the RH side of 
each component of the column vector on the LH represents a different final (outgoing) 
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+ 
en per? p? Si Si2 S13 Sij €1 1;€ 2 2 
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é we , 
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+ $3, S32 S33 $3) + 

él p? E23 ee 2 r€ 2 p? 
Spr Spr Spy vv Spo ee 
t 
€.2 p' »€ 2 p fi ©, p! E12 9? 
et iat Sits Bee Mi Seite, cat etek Skah cas Cn p72 4! 
‘i i 
Final Eigen States S Matrix Initial Eigen States 


eigenstate. 


The Probability of a certain transition occurring is equal to the square of the absolute 
value of each component of the S Matrix linking an initial and final state. For instance, 
the chance of the first initial state, consisting of an electron and positron with specific 
spins and momenta, interacting to form the second end state, consisting of an electron 


and positron with somewhat different specific spins and momenta, is 
St, So1 = |Sai|? = probability of 1st eigenstate transitioning to 2nd. 


As a result, each component of the S Matrix, S;;, represents a transition amplitude 
for a specific response (scattering event, transition, interaction) between particles in 
certain eigenstates, as defined in section 2.1. From here on out, our job is to figure out 


how to calculate Sy; for any given interaction. 


Keep in mind that the S Matrix is huge, almost infinite, because it includes not only all 
possible types of interactions, but also all possible momenta magnitudes and directions 
for each type (all possible eigenstates), and many (most) of the components of the above 
matrix is equal zero, because many interactions such as (e, y > e,e*) are impossible 
because charge is not conserved.The probability that the particles in the 3°¢ (usually 
multiparticle) final eigenstate will be measured to be unchanged from the initial ones 


is represented by |S33|”. 
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2.3.4.2 Physical Space Visualisation 


Lets visualise physically what is happening for this ,Consider an example of an initial 
e~ and e* annihilating one another to become a virtual photon, which in turn mutates 
into a final ~~ and y*. Each possible outgoing state could have different momenta for 
the two particles (with same total momentum). Each such state would be a differ- 
ent possible (multiparticle) final eigenstate. Each would, in general, have a different 


probability of being measured from the other final eigenstates. Of course, we could 


also have final states of fermion pairs e~e* or T~7*, each with many different possible 
momenta states. So, 

|i=1) |f=1) | f=2) 
i i 








i. 
- + - + - + - + 
€ ppl? eben) => Soper => S11 Cpt Craps) +S €,1 p3? ef a) eee pat Hope Hi ps ) +... 








electron-positron final states muon-antimuon final state 


Hence Soper converts the initial eigenstate into a final general state sum of eigenstates 
where the square of the absolute magnitude of the coefficient of each final eigenstate, 


S';, equals the probability of measuring that final eigenstate. 


2.3.4.3 Visualisation in Fock Space 


In Fock space every eigenstate can be visualized as a separate axis in an infinite (ab- 
stract) dimensional space, the Soper can be visualized as a sort of abstract ”rotation” 
in that space. The initial state vector |i) is ”rotated” by the Soper into a new vector 
with components along the eigenstate basis axes. In Fig. 2-2, we illustrate this for a 


final general state composed of only two component eigenstates. 


This should seem reasonable, as Soper has the form of a function e~® where @ is the 
angle a complex quantity (a phasor) in complex space is rotated through. Also, , we 
see that the norm (abstract ”length” = 1 ) of the state vector is unchanged under the 
operation Soper , 1-e., we have a different final vector but it has the same magnitude. 


This is what happens in rotation. 


By close analogy, a complex number can be represent as z = Re’® with length R. 
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Operating on z with e~” is e~’z = eRe’? = Re~**"?, which has the same length 
R, but is rotated through 6. 


\f=3> |F>= Sa1|f=2> + S31|f=3> 


Fock space 
(most dimensions suppressed) 
Soper \f=2> s 
Assumption: only two 
final eigen states 
li=1> = |f=1> 


Figure 2.2: The action of the S operator Visualised in Fock space 


2.3.5 What is H; in the I.P. of QFT? 


We've seen that all of the elements in the perturbative expansion of the S-Matrix 
include time ordered expressions for the interaction Hamiltonian H;. The interaction 
Hamiltonian Hy, is different for various interaction types, i.e., different kinds of forces. 
Electromagnetic interactions are one thing, weak interactions are another, and strong 
interactions are yet a third. Our current aim is to find it in the most basic instances so 
that we may utilise this characteristics of the interacting hamiltonian to establish the 
unitarity of the S-Matrix later in this chapter’s second last part. Recall that there is 
a general rule that the interaction Lagrangian (or wave equation, or Hamiltonian) for 
states in the S.P. of Relativistic QM is the same as that for fields in the I.P. of QFT. 
For free fields, the H.P. and I.P. are the same, so this is consistent with our prior work. 
Using this rule for quantum electrodynamics or QED) for an example , we can find the 


e/m interaction Hamiltonian QED"! via 
Hy = —Ly = —ep(2)y*A,(2)b(2) QED™ = [Ards 


and for QED, this can be used in the perturbative expansion of the S-Matrix(Dyson 
Series). In the I.P., fields (operators) behave like free fields, as we have seen. Given 
this, we take the fields w, 7, and A” as free fields (even though they are part of the 
interaction Hamiltonian). We leave further exploration for proving the unitarity of the 


perturbation form in the second last section of this chapter. 
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We will find that, similar to 3D volume, considering the integration time for the oper- 
ation of Soper to be very large, effectively infinite, will be advantageous in developing 
the mathematics of QFT. Fortunately, again as we will eventually see, the results using 


this assumption can be used to accurately predict real world experimental outcomes. 


2.3.6 In States and Out States 


There are several interactions that may occur in particle theory. A large number of 
different incoming particles can interact with one other. There are numerous final (and 


frequently distinct) sets of departing particles for each set of entering particles. 


A state made up of multiple non-interacting particles can be thought of as a direct 
product of one-particle states under the inhomogeneaus Lorentz group. We use their 
four-momenta p”, spin z-component (or, for massless particles, helicity) 7, and, because 
we may now be dealing with more than one species of particle, an additional discrete 
label n : for the particle type, which includes a specification of its mass, spin, charge, 


and so on, to label the one-particle states. The general rule of transformation is: 








U(A, a) Vp, 01 n13p2,02,na3" = e(—iau(pi'tpe+~ )) 
Api)° (Apa) ' | 
: ! 2 ie Do DBZ, (W (A, pi) DE, (W (A, p2)) (2.23) 
1/2 lay 
O05" 


x Vapi ,o1 najApa,o},noy 


where W(A,p) is the Wigner rotation , and DY) (W) are the conventional (27 + 1)- 
dimensional unitary matrices representing the three dimensional rotation group. (This 
is for massive particles; for any mass- less particle, the matrix D® (WIA, p)) is replaced 
with 6,/, exp(ia@(A, p)) where @ is the infinitesimal transformation group parameter. 


The states are normalized as 


7 / < / / / . . 
(Woh St ail ins Wcunaseunnes) 


ae aint + permutations (2.24) 
=0 (pi a P1) Oo1.019n' m9 (Py ca P2) 004,02 Onk ne ot 


with the term ’ + permutations’ included to take account of the possibility that it is 
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some permutation of the particle types n{,n5,--- that are of the same species of the 
particle types n1,N2,--- (its sign is —1 if this permutation includes an odd permutation 


of half-integer spin particles, and otherwise +1. .) 


We often use an abbreviated notation, letting one Greek letter, say a, stand for the 
whole collection py, 01,71; p2,02,2;--- In this notation, the previous equation is writ- 
ten simply as 

(Wy, Va) = 6 (a —a) (2:25) 


with 6 (a’ — a) standing for the sum of products of delta functions and Kronecker deltas 


appearing on the right-hand side of Eq. (2.24). Also, in summing over states, we write 


prone S° [end'n 


N1O1NQ0Q°* 


In particular, the completeness relation for states normalized as in Eq. (2.25) reads 


U= [eve (U,V) 


The transformation rule (2.23) is only possible for particles that for one reason or 
another are not interacting. Setting A4¥ = 6”, and a” = (0,0,0,7), for which U(A, a) = 
exp(tHr), Eq. (2.23) requires among other things that UV, be an energy eigenstate 


with an energy equal to the sum of the one-particle energies 
Byatt 


and with no interaction terms, terms that would involve more than one particle at a 


time. 


On the other hand, the transformation rule (2.23) does apply in scattering processes at 
times t — -+too. As explained in the introduction in the typical scattering experiment 
we start with particles at time t + —oo so far apart that they are not yet interacting, 
and end with particles at t + +00 so far apart that they have ceased interacting. We 


therefore have not one but two sets of states that transform as in Eq. (2.23): the ‘in’ 
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and ‘out’ states V,* and Y,~ will be found to contain the particles described by the 


label a if observations are made at t + —oo or t + +o, respectively. 


Note how this definition is framed. To maintain manifest Lorentz invariance, in the 
formalism we are using here, state-vectors do not change with time - a state-vector V 
describes the whole spacetime history of a system of particles. (This is known as the 


Heisenberg picture, in distinction with the Schrodinger picture, where the operators 





are constant and the states change with time.) Thus we do not say that V,.~ are the 


limits at t > + of a time-dependent state-vector W(t). 


However, implicit in the definition of the states is a choice of the inertial frame from 
which the observer views the system; different observers see equivalent state-vectors, 
but not the same state-vector. In particular, suppose that a standard observer @ sets 
his or her clock so that t = 0 is at some time during the collision process, while some 
other observer O’ at rest with respect to the first uses a clock set so that t/ = 0 is at a 
time t = T; that is, the two observers’ time coordinates are related by t/ = t —t. Then 
if O sees the system to be in a state UV, O’ will see the system in a state U(1,-—T)V = 
exp(—iHr)W. Thus the appearance of the state long before or long after the collision 
(in whatever basis is used by © ) is found by applying a time-translation operator 
exp(—iH7r) with T + —oo or T > +00, respectively. Of course, if the state is really 
an energy eigenstate, then it cannot be localized in time - the operator exp(—iHr) 
yields an inconsequential phase factor exp (—iE£,T). Therefore, we must consider wave- 
packets, superpositions [ dag(a)W, of states, with an amplitude g(a) that is non-zero 
and smoothly varying over some finite range AF of energies. The ’in’ and ’out’ states 


are defined so that the superposition 





exp(-iHr) [ da.g(a) vt — | dac® g(a) Yor 


has the appearance of a corresponding superposition of free-particle states for T< 
—1/AE or rT > +1/AE, respectively 


To make this concrete, suppose we can divide the time-translation generator H into 


two terms, a free-particle Hamiltonian Ho and an interaction V 


H=H)+V 
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in such a way that Ho has eigenstates ®, that have the same appearance as the eigen- 


states Ut and WV; of the complete Hamiltonian 


Hp ®, = EaPa 
(Oy, ®,) = 6(a’ — a) 


Note that Ho is assumed here to have the same spectrum as the full Hamiltonian H. 
This requires that the masses appearing in Hp be the physical masses that are actually 
measured, which are not necessarily the same as the ‘bare’ mass terms appearing in H; 
the difference if there is any must be included in the interaction V not Ho. Also, any 
relevant bound states in the spectrum of H should be introduced into Ho as if they 


were elementary particles.” 


The ‘in’ and ‘out’ states can now be defined as eigenstates of H, not Ho, 








He = E,W 


a a 


Where W;, is the In state and V+ represent the Out State. They satisfy the condition 





[ dave" g(a) vs > [ dee" 9(a)¥o (2.26) 


for T + —oo or T — +00, respectively. Eq.(2.26) can be rewritten as the requirement 
that: 





exp(-iHr) f dag(a)¥s ++ exp (itr) f drg(a)®. 


for T — —oo or T — +00, respectively. This is sometimes rewritten as a formula for 


the ’in’ and ’out’ states: 





U* = (+00), (2.27) 


a 


where 
Q(7) = exp(+iH7) exp (—iHorT) 


However, it should be kept in mind that Q(-oo) in Eq. (2.27) gives meaningful results 


only when acting on a smooth superposition of energy eigenstates. 


One immediate consequence of the definition (2.26) is that the ’in’ and ’out’ states 


are normalized just like the free-particle states. To see this, note that since the left- 
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hand side of Eq. (2.26) is obtained by letting the unitary operator exp (—iHT) act 
on a time-independent state, its norm is independent of time, and therefore equals the 


norm of its limit for T > oo, i.e., the norm of the right-hand side of Eq. (2.26): 














| de-dB exp (i (Ea — Ep) 7) ola)o"(8) (U5, 2) 


= [eaas exp (—7 (Ea — Eg) T) g(a)g* (8) (®,, Ba) 


Since this is supposed to be true for all smooth functions g(a), the scalar products 


must be equal 














(U5, U2) = (Ge, Ba.) = 4(8 — a) 


It is useful for some purposes to have an explicit though formal solution of the energy 


eigenvalue equation 








HU = EW 


a a 


satisfying the conditions (2.26). For this purpose, write the above equation as 











(Eq — Ho) Va* =VV.* 


The operator Ey — Hp is not invertible; it annihilates not only the freeparticle state 
®,, but also the continuum of other free-particle states ®g of the same energy. Since 
the ’in’ and ’out’ states become just ®, for V — 0, we tentatively write the formal 


solutions as ®, plus a term proportional to V : 








U,+ = 6,4 (E, — Hy tie) VUt 


a 


or, expanding in a complete set of free-particle states, 
WiEs O, 

Ut = 6, 4 [os a 
EF. — Eg tie (2.28) 

Te = (Bg, Ve) 














with € a positive infinitesimal quantity, inserted to give meaning to the reciprocal of 
E,, — Ho. These are known as the Lippmann-Schwinger Equation, which is in 


Scattering theory equivalent to the Schrodiner equation. 
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It can be shown that Eq. (2.28), with a +ie or —ie in the denominator, satisfies the 


condition 





[dace (a) us > [ dz" (a)®e 


for T —> —oo or T — +00, respectively for an ’in’ or an ’out’ state, respectively. For 


this purpose, consider the superpositions 








Ue (t) = . dae *®nt g(a) Ut (2.29) 


6,(t) = [ dac™G(a)o, 


We want to show that U,*(t) and ©, (t) approach ®,(t) for t + —oo and t > +00, 
respectively. Using Eq. (2.28) in Eq. (2.29) gives 


—iEat 
a at a)T gq + Bg 
We (t 
g | @)+ fda fas te: 


Let us recklessly interchange the order of integration, and consider first the integrals 


@- Ps'g(@) Tae 
Zat= /d . 
B i a ee 











For t — —oo, we can close the contour of integration for the energy variable Ey in 
the upper half-plane with a large semi-circle, with the contribution from this semi- 
circle killed by the factor exp (—iE,t), which is exponentially small for t — —oco and 
Im £, > 0. The integral is then given by a sum over the singularities of the integral 
in the upper halfplane. The functions g(a) and T,+ may, in general, be expected to 
have some singularities at values of E, with finite positive imaginary parts, but just 
as for the large semi-circle, their contribution is exponentially damped for t — —oo. 
(Specifically, —t must be much greater than both the time-uncertainty in the wave- 
packet g(a) and the duration of the collision, which respectively govern the location of 
the singularities of g(a) and Tg,+ in the complex E, plane.) This leaves the singularity 
in (Eq — Eg + ie), which is in the upper half-plane for .4;— but not .%3+ We conclude 
then that Fe vanishes for t —- —oo. In the same way, for t — +oo we must close 


the contour of integration in the lower half-plane, and so .%g~ vanishes in this limit. 





We conclude that V5(t) approaches ®,(t) for ¢ + Foo, in agreement with the defining 
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condition 





[dace 9(a) 3 > [ dz" (a)®e 


2.3.7 S-Matrix defination from In-Out States 


An experimentalist generally prepares a state to have a definite particle content at 
t — —oo, and then measures what this state looks like at t ~ +oo. If the state is 
prepared to have a particle content a for t + —oo, then it is the ’in’ state U7, and if 
it is found to have the particle content 8 at t + +00, then it is the out’ state V,—. 
The probability amplitude for the transition a — ( is thus the scalar product 


Spa = (Ug, U2) (2.30) 


This array of complex amplitudes is known as the S -matrix. If there were no interac- 
tions then ’in’ and ’out’ states would be the same, and then $3, would be just d(a—/). 
The rate for a reaction a — (3 is thus proportional to [Sq — 6(a — )|’. We shall see in 


detail in subsection 2.3.10 what Sq has to do with measured rates and cross-sections. 


Perhaps it should be stressed that ‘in’ and ’out’ states do not inhabit two different 
Hilbert spaces. They differ only in how they are labelled: by their appearance either 
at t + —oo ort > +00. Any ”in’ state can be expanded as a sum of out’ states, with 
expansion coefficients given by the S -matrix (2.30). Since Sig, is the matrix connecting 
two sets of orthonormal states, it must be unitary. We Il Prove its Unitarity in the last 


section and also discuss some of its implications. 


2.3.8 Analytical form of S -Matrix(Born Approximation) 


It is often convenient instead of dealing with the S -matrix to work with an operator S, 
defined to have matrix elements between free-particle states equal to the corresponding 


elements of the S -matrix: 
(Pz, S®,) = Be 


Recall the explicit though highly formal expression for the ’in’ and ’out’ states i.e 





We, = oo) Dy, 


a 
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yields a formula for the S -operator as : 
S = 2(00)'N(—00) = U(+00, —00) 
where 
U (7,70) = Q(7)'Q (79) = exp (iHoT) exp (—iH (7 — 7)) exp (—i HoT) 


This will be used in the next section to examine the Lorentz invariance of the S -matrix, 
a The methods of the previous section can be used to derive a useful alternative formula 


for the S -matrix. Let’s return to Equation 


—iEat 
L at a)T gq + Bg 
We (t 
g | @)+ fda fas ae: 


for the ’in’ state UT, but this time take t — +oo. We must now close the contour of 








integration for EF, in the lower hatf- E, -plane, and although as before the singularities 
in T3,* and g(a) make no contribution for t + +-oo, we now do pick up a contribution 
from the singular factor (E, — Eg + ie)”. The contour runs from E, = —oo to E, = 
+oo, and then back to E, = —oco on a large semi-circle in the lower half-plane, so 
it circles the singularity in a clockwise direction. By the method of residues, this 
contribution to the integral over E, is given by the value of the integrand at Ey = 
Eg — te, times a factor —2iz. That is, in the limit « — 0+, for t + +00 the integral 


over @ in equation has the asymptotic behavior 
Int + —Line ‘Pet / did (Eq — Eg) g(a)T3a* 
and hence, for t + +00, 
tt) [epee a = din fdas (EF, — Eg) g(a)Tga* 


But expanding the superposition 








v(t) = [ dac™'g(ayus 
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for U,* in a complete set of ’out’ states gives 


We (i [eee tg(a) f a8 V5 S30 


Since Sgq contains a factor 6 (Eg — E,), this may be rewritten 


VF (t) = [ asuzet® f dag(a) Sp 


and, using the defining property (2.26) for ’out’ states, this has the asymptotic behavior 
for t 4 +00 


W(t) [ a3@ae%* [ dag(a)S 


Comparing this with our previous result, we find 


| 4e9(0) Sin = (8) — 2in fxs (Eq — Eg) g(x) TH, 


or in other words 
Se _ 6(8 = a) — 276 (Ea = Ez) Daa 


This suggests a simple approximation for the S -matrix: for a weak interaction V, 


we can neglect the difference between ‘in’ and free-particle states in the expression 








To = (®g, VW) in which case Sgq becomes 
SBa = d(C = a) — 2176 (Eo a Ez) (Pz, V®,) 


This is called Born Approximation.Which is reasonably good for very weak inter- 
actions.For higher order terms we had to turn to the Perturbative expression of the 
S-Matrix. 


2.3.9 Symmetries of S-Matrix 
In this section we will consider both what is meant by the invariance of the S-matrix 


under various symmetries, and what are the conditions on the Hamiltonian that will 


ensure such invariance properties . 
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2.3.9.1 Lorentz Invariance 


For any proper orthochronous Lorentz transformation x — Aw +a, we may define a 


unitary operator U(A,a) by specifying that it acts as follows 
UA, A) Vp, 01 n1sp2,02,nos~ = exp (—ia, (pt A: Pb Bi )) 


Xx [om (Apa) ~ Do (Ww p1)) De (W (A, p2)) 2 








PiPo°** 


O05 
x W api .o4 na sAp2,o4,nos 
on either the ’in’ or the ’out’ states. When we say that a theory is Lorentz-invariant, 


we mean that the same operator U(A,a) acts as above on both ’in’ and ’out’ states. 


Since the operator U(A, a) is unitary, we may write 
Sea = (Wz, 8) = (U(A, a) Up — U(A, a) 82) 


so using (2.23), we obtain the Lorentz invariance (actually, covariance) property of the 


S -matrix: for arbitrary Lorentz transformations A“, and translations a“,as 


OS pl Daas coh: ipl ‘ Ried 
Spf of nt 3P'21FQ)MQ5°"* 7P15,01,N15P2,02,N2) 


= exp (ia, (pi + ph +--+ — pit — pf! —---)) 
x / (Api)° (Apa)! «- (Api)? (Ape)? = 








PIP pepy ++: 
x 2 DE (W(A, p1)) D2, (W (A, p2)) = 


0102" 
/ 
. * J * 
(31) (1 (sunt) ) DD, (w+ (4.06))~ 
x y Si 
a5 
x Spt at ni :Aph,a$ Nb ,Ap1,01 ni; Ap2,02,n23-" 
(Primes are used to distinguish final from initial particles; bars are used to distinguish 
summation variables.) In particular, since the left-hand side is independent of a”, so 


must be the right-hand side, and so the S -matrix vanishes unless the four-momentum 


is conserved. We can therefore write the part of the S -matrix that represents actual 
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interactions among the particles in the form: 
Spa — 6(8 — a) = —27iMg,.5° (pp — Pr) 


(Note: the amplitude Mg, itself contains terms that involve further delta function 


factors. ) 


Above equation should be regarded as a definition of what we mean by the Lorentz 
invariance of the S -matrix, rather than a theorem, because it is only for certain special 
choices of Hamiltonian that there exists a unitary operator that acts as in (2.23) on 


both ’in’ and ’out’ states. 


2.3.9.2 Internal Symmetries 


There are various symmetries, like the symmetry in nuclear physics under interchange 
of neutrons and protons, or the ’charge-conjugation’ symmetry between particles and 
antiparticles, that have nothing directly to do with Lorentz invariance, and further 
appear the same in all inertial frames. Such a symmetry transformation JT’ acts on 
the Hilbert space of physical states as a unitary operator U(T), that induces linear 


transformations on the indices labelling particle species as 


UE) Upawinipseraye = oy Bran. (T) Fraing (T) i. 
AN (2.31) 


x Vproiai spoon: ‘S? 


In accordance with the general property of Transformation Group the U(T) must satisfy 


the group multiplication rule 


where TT is the transformation obtained by first performing the transformation T, 
then some other transformation T. Acting on Eq. (2.31) with U(T), we see that the 


matrices J satisfy the same rule 
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Finally, taking the scalar product of the states obtained by acting with U(T) on one 
out’ state and one ’in’ state shows that Y commutes with the S -matrix, in the sense 
that 


* = on a 
. : ; ps int Ném (TL Prin (LT )DranglP) x Spiot.N 137505,.N3;-- ,p101.N1;p202.No;-- 
Ni No. NING. 
= = Spi ot nh phosnhso1ni; }p202N2,"° . (2.32) 


Again, this is a definition of what we mean by a theory being invariant under the 
internal symmetry 7, because to derive Eq. (2.32) we still need to show that the 
same unitary operator U(T) will induce the transformation (2.31) on both ’in’ and 
‘out’ states. This will be the case if there is an ’unperturbed’ transformation operator 


Uo(T) that induces these transformations on free-particle states, 


Uo(T) Pprornripeoana: =) De wll )Prigng (LT ) aXe D ecivipaa Nets: 
N, No... 


and that commutes with both the free-particle and interaction parts of the Hamiltonian 


Up (TL) HoUo(T) = Ho 
UT WUy(T) =V 


From either the Lippmann-Schwinger equation (2.27) or from (2.28), we see that the 
operator U,(T) will induce the transformations (2.31) on ’in’ and ’out’ states as well 
as free-particle states, so that we can derive Eq. (2.31) taking U(T) as Uo(T). 


A special case of great physical importance is that of a one-parameter Lie group, where 


T is a function of a single parameter 6, with 


In this case the corresponding Hilbert-space operators must take the form 


U(T(0)) = exp(iQ0) 
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with Q a Hermitian operator. Likewise the matrices 0(T) take the form 
Drin(L(O)) = dn exp (idn9) 


where q,, are a set of real species-dependent numbers. Here Eq. (2.32) simply teils us 


that the qs are conserved: Sz. vanishes unless 
Qn’, + Oni toe: = Gn + Ing T°** 


The classic example of such a conservation law is that of conservation of Charge—(Which 
is an intrinsic property of any particle), its conservation denotes a specific case of con- 


servation of Internal Symmetry. 


2.3.9.3 Parity 


To the extent that the symmetry under the transformation x — —x is really valid, 
there must exist a unitary operator P under which both ’in’ and ’out’ states transform 


as a direct product of single-particle states: 





PMS es risieseunon = Mny"n2 °°" WV Bpio1n1;Pp202n2, 


where 7, is the intrinsic parity of particles of species n, and # reverses the space 
components of p“. (This is for massive particles; the modification for massless particles 


is obvious.) The parity conservation condition for the S -matrix is then: 


Spt on! sphohnhs- proanispaoanas 

a * * 

X Spy ont: Ppsohnhs-,Pproiny;Pr2022} 
Just as for internal symmetries, an operator P satisfying the Lorentz transformation 
of S-Matrix will actually exist if the operator Pp which is defined to act this way on 


free-particle states commutes with V as well as Ho. 
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2.3.9.4 Time Reversal 


It can be Proven that that the time-reversal operator 7 acting on a one-particle 
state Upon gives a state Vy) _~>n» with reversed spin and momentum, times a phase 
¢n(—1)2~?. A multi-particle state transforms as usual as a direct product of one-particle 
states, except that since this is a time-reversal transformation, we expect ‘in’ and ’out’ 


states to be interchanged i,e 








Pe ees: = Canty hag —< Wop (-o1)rasPoa(—oa) ha, (2.38) 


Above expression can be expressed in short hand notation as 








TUE = UE 


a 


where 7 indicates a reversal of sign of three-momenta and spins as well as multiplica- 


tion by the phase factors shown in Eq. (2.33). Because T is antiunitary, we have 
(W5, U2) = (TVZ,TV;) 
so the time-reversal invariance condition for the S' -matrix is 
Sa = S70,7B 


or in more detail 


Sy opt ueuln iumtbeaaie = Gall! Ga AIS eee =I a1 eres 


x5 gap, (—01)m1;Pp2(—02)24-. PP, (01) Pho) Jy 


Note that in addition to the reversal of momenta and spins, the role of initial and final 
states is interchanged, as would be expected for a symmetry involving the reversal of 


time. 


The S -matrix will satisfy this transformation rule if the operator 7p that induces 


time-reversal transformations on free-particle states 


To®g = Bea 
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commutes not only with the free-particle Hamiltonian (which is automatic) but also 


with the interaction: 
T) Holo = Ho 


T) Vio=V 


2.3.9.5 Charge Conjugation 


As already mentioned, there is an internal symmetry transformation, known as charge- 
conjugation, which interchanges particles and antiparticles. Formally, this entails the 


existence of a unitary operator C, whose effect on multi-particle states is: 








Cut = Cribs rare |S oe (2.34) 


P101N1 ;p202N2;"- : pion’ ;spooan$,-- 
where n° is the antiparticle of particle type n, and €, is yet another phase. If this is 
true for both ’in’ and ’out’ states then the S -matrix satisfies the invariance conditions 
a? * * 
Spi on) sphobnys- P1O1N1;P2T2N2Q;""° (Siena sO Eni Ene fare pba ieltph oune pe he psaaness 


As with other internal symmetries, the S -matrix will satisfy this condition if the 
operator Cy that is defined to act as stated in Eq. (2.34) on freeparticle states commutes 


with the interaction V as well as Ho; in this case, we take C' = Co 


2.3.10 S-Matrix as a tool to find Transition/Decay rates and 


cross sections 


The S -matrix Sq is the probability amplitude for the transition a + £6, but what 
does this have to do with the transition rates and cross-sections measured by experi- 


mentalists? In particular, 
Spa — 6(8 — a) = —27iMgad* (pg — Pe) 


shows that Sq has a factor 6* (pg — p,), which ensures the conservation of the total 
energy and momentum, so what are we to make of the factor [54 (pg —pa)|’ in the 


transition probability |S Pel The proper way to approach these problems is by studying 
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the way that experiments are actually done, using wave packets to represent particles 
localized far from each other before a collision, and then following the time-history of 
these superpositions of multi-particle states. In what follows we will instead give a quick 
and easy derivation of the main results, actually more a mnemonic than a derivation, 
with the excuse that no interesting open problems in physics hinge on getting the fine 


points right regarding these matters. 


We consider our whole system of physical particles to be enclosed in a large box with 
a macroscopic volume V. For instance, we can take this box as a cube, but with points 
on opposite sides identified, so that the single-valuedness of the spatial wave function 
requires the momenta to be quantized 


27 


p= Fy. (n1, N2, Ng) 


where the n,; are integers, and L? = V. Then all three-dimensional deltafunctions 


become 
/ 


df (p’ — 





oe V 
3,,,0(p—p’)x __ 
[ie re = Bg 


where dp» is an ordinary Kronecker delta symbol, equal to one if the subscripts are 


p)= Ons 


equal and zero otherwise. The normalization condition (2.24) thus implies that the 
states we have been using have scalar products in a box which are not just sums of 
products of Kronecker deltas, but also contain a factor [V/ (27)3]”, where N is the 
number of particles in the state. To calculate transition probabilities we should use 


states of unit norm, so let us introduce states normalized approximately for our box 


No/2 


Wem = [(2n)2/V]" We (2.35) 


with norm 
(whe eo ) = bBo 


where 0gq is a product of Kronecker deltas, one for each three-momentum, spin, and 
species label, plus terms with particles permuted. Correspondingly, the S-matrix may 


be written 


Spa = [V/(2n)3] NetNe)/? gBox. 


where $3 is calculated using the states in equation (2.35). Of course, if we just leave 


our particles in the box forever, then every possible transition will occur again and 
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again. To calculate a meaningful transition probability we also have to put our system 
in a ‘time box’. We suppose that the interaction is turned on for only a time 7’. One 


immediate consequence is that the energy-conservation delta function is replaced with 


1 pte 
br (Ba — Ep) = 5 / ane i iar 
1 
The probability that a multi-particle system, which is in a state a before the interaction 


is turned on, is found in a state 6 after the interaction is turned off, is 
Ox NatN, 2 
P(a + B) = [SB]? = [(2n)9/V]™*) |S 


This is the probability for a transition into one specific box state @. The number of 
one-particle box states in a momentum-space volume d?p is Vd*p/(27)?, because this 
is the number of triplets of integers n1,n2,n3 for which the momentum 

27 


a (n1, no, n3) 


ae 


lies in the momentum-space volume d’p around p. We shall define the final-state 
interval d3 as a product of d*p for each final particle, so the total number of states in 
this range is 


dNp = [V/(20)*] ap 


Hence the total probability for the system to wind up in a range d( of final states is 


dP(a > 8) = P(a > B)dNg = [(20)3/V}"* |S, |? dB (2.36) 


We will restrict our attention throughout this section to final states 6 that are not 
only different (however slightly) from the initial state a, but that also satisfy the more 
stringent condition, that no subset of the particles in the state (other than the whole 
state itself) has precisely the same four-momentum as some corresponding subset of 
the particles in the state a ,this means that we are considering only the connected part 
of the S-matrix.) For such states, we may define a delta function-free matrix element 
Mee: 
Spa = —2indy, (Ps — Px) Or (Es — Ea) Mpa 


Our introduction of the box allows us to interpret the squares of delta functions in 
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|Sgal? for 6 4 a as 


[5% (Ps — Pa)] = 88 (Ps — Pa) 53. (0) = 53 (Ps — Pa) V/(2m) 
[ir (Eg — Ea)|? = br (Eg — Ea) 6r(0) = 67 (Eg — Eq) T/20 


so Eq. (2.36) gives a differential transition probability 


dP(a + B) =(2n)? [(20)8/V]"** (T/2m) |Moal” 


x Oy (Pg i Pa) or (Eg = Eu) dB 


If we let V and T' be very large, the delta function product here may be interpreted as 
an ordinary four-dimensional delta function 6* (pg — pa). In this limit, the transition 
probability is simply proportional to the time 7’ during which the interaction is acting, 


with a coefficient that may be interpreted as a differential transition rate: 


dl'(a > 6B) =dP(a > B)/T 


3Na—2171-N, 2 <A (2.37) 
= (20) VO" |Mgal O° (Pps — Pa) a8 


where now 
Spa = —27id" (pg — Pr) Moa (2.38) 


This is the master formula which is used to interpret calculations of S— matrix elements 
in terms of predictions for actual experiments. We will come back to the interpreta- 
tion of the factor 64 (pa — pg) dG later in this section. There are two cases of special 


importance: 
No =1: 


Here the volume V cancels in Eq. (2.37), which gives the transition rate for a single- 


particle state a to decay into a general multi-particle state 6 as 
dV (a > B) = 20 |Mgo|” 5* (ps — Pa) dB 


Of course, this makes sense only if the time 7’ during which the interaction acts is 
much less than the mean lifetime T, of the particle a, so we cannot pass to the limit 
T > oo in by (Ey — Eg). There is an unremovable width AE ~ 1/T 1/7, in this delta 
function, so above equation is only useful if the total decay rate 1/7, is much less than 


any of the characteristic energies of the process. 
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Ne = 


Here the rate (2.37) is proportional to 1/V, or in other words, to the density of either 
particle at the position of the other one. Experimentalists generally report not the 
transition rate per density, but the rate per flux, also known as the cross-section. The 
flux of either particle at the position of the other particle is defined as the product of 
the density 1/V and the relative velocity ua : 


DO, = tia Vv 


(A general definition of uy is given below; for the moment we will content ourselves 
with specifying that if either particle is at rest then u, is defined as the velocity of the 


other.) Thus the differential cross-section is 
do(a > B) = dl (a > B)/®, = (27)*uz* |Mgal? 54 (pg — pa) dB 


Even though the cases Na = 1 and Ng = 2 are the most important, transition rates 
for Na => 3 are all measurable in principle, and some of them are very important in 
chemistry, astrophysics, etc. (For instance, in one of the main reactions that release 
energy in the sun, two protons and an electron turn into a deuteron and a neutrino.) 
subsection 2.3.11.2 presents an application of the master transition rate formula (2.37) 


for general numbers N, of initial particles. 
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2.3.11 Unitarity Of S-Matrix in Canonical QFT’s and Impli- 


cations 


In elementary particle physics, the S -matrix plays a fundamental role for describing 
scattering processes. In what follows, we want to describe the basic ideas of this 
fundamental concept in terms of a complex finite-dimensional Hilbert space X which 


has the dimension N = 1,2,... Suppose that the two systems 


P1,---PN and %1,...,Wn 


form an orthonormal basis of the Hilbert space X. By the Dirac calculus, we get the 


following key relation: 


N 
So Wi | er) (er | Ue) = | Ue) = Oe, R=... 
r=1 


S -matrix. By definition, the transition amplitude from the state ~, to the state y; 
reads as 
Sin = (9%; | Ye), GHA 


The idea is to regard this system (s,,) of transition amplitudes as the elements of a 
complex (N x N) -matrix called the S -matrix as we have seen in subsection 2.3.4. 
From above equation we obtain immediately the following crucial property that The S 
-matrix (s;,) must unitary. Explicitly, 

N 

So shi Sirk =Ojn, j,k=1,...,N 


p=, 


Transition probabilities. Letting 7 = k, we get 


N 
» bal Sle BS lal 
r=, 


This allows us the following interpretation. For given fixed state w, with k = 1,...,N, 


the number 
sal = Oe (Mig FHS Ley 
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is the transition probability from the state wv, to the state y,. We now want to show 


that precisely unitary operators are related to systems of transition probabilities. 


From transition probabilities to unitary operators: Let y1,...,gy and ¥j,...,vUN 


be chosen as above. Define a linear operator S': X + X by setting 
Spr = We; hee ete sad 
Then, the operator S is unitary, that is, 


(Sp| Sv) =(~|v) forall g,pex 
Note first that 
(Sip; | Ser) = (by | ti) = bf, b= 1,...,N 


Thus, it follows from y = S7™_, (ye | ¢) ye and = 3%, (yr | Y) yr that (Sy | Sx) is 
equal to 


N N 
Swe |v)! (er |b) (Sex | Ser) = 3° | Ge) (oe | H) = (Y | Y) 


kl=1 k=1 


From unitary operators to transition probabilities: Conversely, suppose that 
we are given a linear unitary operator S: X > X. Let y1,...,yy~ be an orthonormal 
basis of X. Finally, set 

Wp tna. ‘US ee 


Then, the system w,...,#y is also an orthonormal basis. In fact, 


(We | Vi) = (Ser | Sei) = (er | Yi) = Oye 


The matrix elements of the operator S are given by 
Si = (yO OR)s eS Les lh 


Hence sj, = (y; | Uz), Showing us that sj, represents a transition probability, as in 


above expression. 


Physical interpretation. Suppose that the following situation is met. - The given 


quantum system is in one of the states y1,..., yy at the initial time ¢;, . - The initial 
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state y, passes over to the state Sy, at the final time tout . 


Then, the real number 
Issel = (v5 | See)? 


represents the transition probability of the quantum system from the initial state yp, 
to the final state y; during the time interval [tin , tout |. It is crucial that the operator 
S is unitary. Otherwise, an interpretation in terms of transition probabilities is not 


possible. 


Typically, the original S -matrix constructed in the Standard Model in particle physics 
is not unitary. However, this shortcoming can be fixed by introducing ghosts and 


antighosts. 


2.3.11.1 Proof of Unitarity of S-Matrix in Canonical formalism 


1.Unitarity of the Analaytical form : Recall that the completeness of in-out states 


relation states 


w= f da|¥o) (a | ¥) 


and the requirement of in-out states is that they should be a energy eigenstates satis- 


fying the Schrodinger equation 


with 
Ey =pi+petpe+.... 


where p{, p5..are energies of single particles and E is the sum of all single particle 


energies. The Hamiltonian can be split into two parts as 
H=H)+H, 
Where Ho is the free Hamiltonian satisfied by 


Ho®q = Ea Bq (2.39) 
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and ®, has the same appearance as in-out states as discussed in section 2.3.6 and they 


are orthonormal as 


(®, | ®g) = d(a — B) 


or 


(Oy | Pa) = 4 (a’ <0) 


The in-out states are defined as eigenstates of total Hamiltonian H 
H=H +H, 


i,e 








Hwv= = Ev 


a a 


and they satisfy the equation 





[doc g(ayus > [doc G(a)o, 


T —> —0o,T —> +00( free states in long past and long future) and 





cit) F daglay vs — eit) F dagla)®, 


for T 4 —00o,T + +00 

















Hence 
WS = Q[boo) ®, 
where 
Q(+00) = e(ttHr) 6(-iHor) 
Also 
(U5, Uz) = (Sz, ®,) = d(8 = @) 
since 








HU = EW 


a a 








=> (Hy + H) We = E,W 


a 











—> Hyvt — E,Wt = —H,vt 


Qa 
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Or 











(Eq — Ho) UZ = HV 


and formal solution of this will be 








wt = 6, + (Ey — Hy tie) | Hy wt 


In Dirac notation 








|W*) =|.) + (By — Ho + ie) Hy |V2) 


a 


Inserting the completeness relation 


i: d8|®,)(®,| =1 








— |v) =|8,) + i d8|®3)(®5| (Eq — Ho + ie)! Hy |W) 








1 £4 
==. jb.) + f dB |os) (E. — Eg + ie) (®g |H;| V2) 








- 
= [6+ f ab, pa__1,) 








E., — Ep + ie) 
where 
Taq = (|r| Va) 
In-state particle content = state a; at t = —oo and out state particle content = state 


B at ; t = +00 Now the S-Matrix is denoted by 
Spa = (U5 | 2) 


NOW Wwe Call See 


| 40SiySsu— f a30vs |W" (a | ME) 
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= / dB(WE)W5) (U5 | WE) = (WE i dB |W5) (Vahey = (wtpwe) 


since 














(U5 1a) = 6(6 — a) 


hence => (WF Wz) = 6(y — a) and finally we get 


/ dB5%,,S¢a = 6(y — a) > StS =1 


Similarly if 
[48519835 =i(y-a) 3 SSt=1 


Proof: According to Born approximation 


Spa = 5(B — a) — 20d (Eq — Ep) Tae, 
where 
6 (Eq — Es) = Energy conserving delta function 
6(6 — a) =Trivial amplitude factor 
and T;, = (®,|H7|V¥2) 
Actually for Born approximation 


SBa ~ 6(6 = a) — 2776 (Eo > Ez) (®,|H7|®.) 


ie |U*) = |®,)(approximated as free particle state). 


Now consider the Most general Matrix element given by 








(U5 |Ai| We) 


Recall we had 








|\W=) = |®,) + (EZ, — Ho + ie) ' A; | VF) 
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then 








(Wg*| = (®s| + (WF (Es — Ho £ ie) * Hj 


since 














(Us* |Ar| Va) = (Ue |Ai| We) 








hence putting the expressions for |W*) and (W,*| we get 


Qa 














(5 [Hy ([Ba) + (Ba — Ho ie)" Hy |WE)) = ((-+ (WEI (Np Hokie Hi) M2) 








=> (Vg + |Hy| ®,) + (UF | Ay (Ey — Ho + ie)" H,| VE) 
= (8g |Hz| Vat) + (5 |(Eg — Ho £ ie) ' H3| HV.+) Nowwehave(2.40) 

















(UF |H,H,| VE) 


(Ey — Hy ie) © 


(Wo+Hy*H,U+) 
(E3 — Hy + ie) 





Summing over complete Set of |®,) states we get 


flo. =1 














Up Ar J dy] ©) (Py LAr Ya) | (Ue Hi S dy | Gy) (@y | Hr Va) 
| (Eg — E, + ie) 

















2K T _ 
ne (E, — Ey + ie) 














— f dy (WE | Hi®,) (®, | HyWS) fy (WEHFS,) (®, | HW 
(E. — E, £ ic) (E3 — E, £ie) 


os 
































or 


























Ty =- fay TB 3 


x ([Ea — Ey tie] — [Eg — E, #ie)™*). 


(2.41) 


To prove the orthonormality of the ’in’ and ’out’ states, divide by Ey — Eg + 2ie. This 
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gives 











Ty ; Tie 
Eg—Ey+2ie}] © Ey — Eg + ie 


= - fey Te Tyat 
E,—Eqytie} Ey — Eg tie 


The 2es in the denominators on the left-hand side can be replaced with es, since the 








only important thing is that these are positive infinitesimals. We see then that 6(3 — 
a) + T3,/ (Es — Ey £ ie) is unitary. With Lippmann Schwinger Equation(2.28), this 


is just the statement that the V> form two orthonormal sets of state-vectors. 








The unitarity of the S-matrix can be proved in a similar fashion by multiplying (2.41) 
with 5 (Eg — E,) instead of (Ey — Eg + 2ie)™*. 


Multiplying both sides by 6 (Eg — Eq) 























T:46 (Ey — Ey) -T#6(Es — Ea) = i dy 
(2.42) 

















Since according to born approximation 

Spa = 5(B — a) — 20id (Eq — Ep) Tae, 
Grouping the terms with similar coefficients we finally get 

| 485348 =é(y-a) > S'S =1 


Hence S -Matrix is always Unitary 
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2.Unitarity of the Perturbative/Dyson Series form : In section 2.3.1 we have 
derived the perturbative form of S-matrix as Dyson series expansion , which is given 


as 





es 3 i i‘ i di: 7 7 dt T[Hy(ts)...Hi(tn)] (2.43) 


Where H7(t,,) is the Interacting Hamiltonian .We convert this into Hamiltonian Density 


since H;(t,,) is more useful in field theory. 


In terms of Hamiltonian Density we have 


Ay(tn) — [ bent) 


Hence the expansion becomes 





so 3 (-i)" / Pte i * Py T Br (ts)Hilto).--Mr(te)] (2.44) 


n=0 oo —oo 


where 7[H7(t1)H7(t2)...H7(t,)] denotes the Time ordered product of the Interacting 


Hamiltonian described at different point in time and here t; > ty > t3... > th. 


Now we have already described that for a theory to be renormalizable , the interacting 
part cannot be more that ¢* in terms of field variable. Since the typical interacting 


Hamiltonian already has a certain ¢? term, hence in the integral 


ie d*tnT [Hr (t,)H, (t2)| 


(oe) 


only 2 terms of H;(t,).H7(t2) suffices. Since the value of coupling constant which is 


denoted as \ and for a typical H,(t) it becomes 
H,(t)=A¢* or Ad* 


Then, 
H1(t1).Hr(te) ~ Ne 


which is smaller than \. 


Similarly for three terms 


H7(t1).Hy7(te)Hr(ts) oie 
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which is very small compared to. 


Hence, we will consider upto the 2nd order terms of S-Matrix expansion. We have 








(iY f= * 
Sh = d al - d‘a1... i- d*tnT (H1(ti)H1(t2)| (2.45) 
If we expand the individual order term as 
(oe) 1 too too 

Sy=1-i f Hadley — 5 f / TiMs(tr)Hy(te)|dtayd42o +. (2.46) 

Therefore ss 
S=S+4+ 594594... = 5% 5 (2.47) 

n=0 


We take only upto 2nd order ,hence 
S=8%94 9% 4 50) (2.48) 


also, 
Ot = SiO) “te gi@) 4 gi) (2.49) 


For Unitarity STS = 1—= SS then 


StS ~ | StO 4 91) 4 st) Is 49% 4 g@) 


~ St gO 4 st gM 4 gt 2 4 gt) gO (2.50) 
4+ ST) 9 4 gt) G2) 4 gt) 90 4 S12) gM 4 gt) g2) 
Hence,we get 
Sth x14 OSM 4 4 st@g@) (2.51) 


Since we already have 1, which implies that the higher order multiples must vanish. 


Let us try to prove this by taking the most general last term St?)S@) We have seen 
that how S-Matrix creates a asymptotic state @; in the long part t = —oo to another 
asymptotically free state @¢ in the long past. That is matrix elements of 5 can be 


written as 


Si = (@f|S|:) 
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However there are infinitely many states between ¢; and @». Now for unitarity of the 


S-Matrix, we have seen that in general for 2 intermediate steps say, 
SS seein. jeSIN 
r=1 


set j=k then 
T = 
pain a 
r=1 


Let j=i( Initial state) and r=f(Final State) we get 


ye S18 =1 
P=L 


Similarly we have 


Shi = (b5|S|d4) 


and 
St. = (di|S*1d5) 
then we get 
5 57St, =1 
i=1 
Hence ss 
Yo (bs 1S1¢i)(b:1 S167) = 1 (2.52) 
i=1 
=> (6/|5 916i) (Gi S*1¢s) = 1 (2.53) 


i=1 


By using completeness relation we can write 


(oe) 


S516: (ei = 1 


t=1 


so we can write 


(o,|SS"|b7) = 1 (2.54) 
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Therefore we can say that,in the expansion of Dyson series 
Sstw14+ SV st 4 4 s@ st) (2.55) 


all other terms must vanish(Except the lst term),if the S-Matrix has to be Unitary. 


For the last term we have 


1 CO CO 
$(2) = -= | ‘i T[Hr(ti)Hr(t2)|d*ad* x. ty > te (2.56) 
1 CO CO 
Si) = -= | is T[Hi(ti)H1(t2)|'d*0,d*x2 t1 > to (2.57) 
= Se as T{H'(t2)H'(t,)|d*x,d* ist 2.58 
ae p(ta)H7(ti)|\d° aid’ ae 1 > te (2.58) 


Let H\(t2) = A(x) and H!(t,) = B(y), we already know that 
T[A(x)B(2)] = O(a — y)A(x) B(y) + O(y — «) By) A(z) 
since ty > t2,s0 y > x. Then only 2nd term survives. 
T[A(x)B(y)] = Oy — 2) By) A(z) 


T [HI (to)! (t1)] = HE (tr) HY (te) 


Then finally 
1 CO CO 


Put these in equation 2.50 and we get 


(¢,|(SOGTO 4 4 SOSA g-) = 1 (2.60) 
1 sis 1 we 
(9,|SOSI +. 4(-5 / H(t) Hr(te)d'aid'x2)(— > / [A (t1)H} (t)|d*x,d*x2)\o) = 1 
(2.61) 


take the last term, that is St) 64) 


- 5 / - / : U1 (ty)d!x1H' (ta)d*x2)|o5) (2.62) 
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ae af Hi (ti)d*2, ye H! (t)d*x2)|o/) (2.63) 


We know that in the interaction picture 


H(t)| W(t) = an Z|WCA 

then it sae a 
ins H(t) | / = ih lop\din| (2.64) 
= 5 [Heater] J inlo;(00)) — Jes(—00)) (2.65) 


|b¢(00)) = |bF) |b¢(—00)) = |bF) 


Hence the above equation becomes zero. that is St@|d-) = 0, so we can write 
SAGA |o-) = 0. 


so in general we can write 
SMS] b,) = 0 (n = 1,2,3,..) 


we can conclude that only the zeroth order product survives,hence 


(b,|SOSt|b7) = 1 (2.66) 
(byl]b7) =1 (2.67) 
(dslo7) = 1 (2.68) 
Hence for SS? = 1 
SS"|b5) = |@F) 
SSt=T 


2.3.11.2 Implication of Unitaritry of S-Matrix 


A. Optical Theorem The unitarity of the S -matrix imposes an interesting and 


useful condition relating the amplitude Maja for forward scattering in an arbitrary 
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multiparticle state a to the total rate for all reactions in that state. Recall that in 
the general case, where the state G may or may not be the same as the state a, the S 


-matrix may be written as 
Spa = 5(8 — a) — 2nid* (ps — Pa) Moo 
The unitarity condition S'S = 1 then gives 
6(7 -— a) = i 485% Sox = 5(-y — 2) — 2r8" (py — Pe) Myo 
+ 2nid* (py — Pe) Mjy + 4” i. dB6* (pg — py) 8° (pe — Pa) Mg, Mga 
Cancelling the term 6(y — x) and a factor 2764 (py — px), we find that for py = pa 
0 = -iMig + iM%, + 2 ) dB5* (pg — Pa) MjyMoa 
This is most useful in the special case a = y, where it reads 


Im Moa = —7 / dB5* (pg — Pa) |Mgal” 
Using the expression of differential transition rate i,e 


dl'(a — 8) =dP(a > B)/T 
= (21)?Ne? V3" No | Meal” 54 (Dg — Pa) dB 


Using the expression for Im M,,. the above equation can be expressed as a formula for 


the total rate for all reactions produced by an initial state a in a volume V 


res fas? 


= (2n)3Ne-2y 1M / 854 (pa — pa) |Meal? 


1 
== (One yt lM 
7 
In particular, where a is a two-particle state, this can be written 


Im Mra = —UgFa/ 167° 
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where Ug is the relative velocity given by 





Ua = V(r po)? = mim3/E\ E» 


where pj, p2 and m1,™mz are the four-momenta and mass of the two particles in the 
initial state a. in state a, and o, is the total cross-section in this state, given by 


integrating 
do(a + B) = dI(a + 8)/®q = (2m)*uz* |Mpol” 5° (pe — Pa) dB 


hence 


a= f dBdo(a > B)/d3 = (2n)*u," f dB [Mpal®5* (pa ~ Pa 


This is usually expressed in a slightly different way, in terms of a scattering amplitude 
f(a — 6). We know that the differential cross-section for two-body scattering in the 
center-of-mass frame is given by 

da (x — £B) (27)*k! Et EE E> 


da hE? IMsal’ 








where k’ and k are the magnitudes of the momenta in the initial and final states. We 


therefore define the scattering amplitude as 


An? /[k/E' ES ELE 
f(a 8) =-— Mea 


so that the differential cross-section is simply 





da(a — £) 


oT = Ifa + 8)? 


In particular, for elastic two-body scattering, we have 


dr? Ey E> 


f(a > B)=- = 


Moa 
now we have seen that if @ is a two-particle state, then 


Im Moa = —UaTo/ 160? 


54 


2.3. S -Matrix and Scattering Theory in Canonical QFT 





if we take pioe 
T 
a> B)= <a 


and put f(a — qa) and take Im f(a > a) we get 


Mea 


An? FE, E 
Im f(a a)= pe iy 
E 
An? Fy E> 
E 


we now use the expression of relative velocity ug given by 


Im f(a a)=- .(—UgOo/ 167°) 





Ua = Vir po)? = mim3/E\ E» 


For Simplicity if we work in the ’center-of-mass’ frame, where the total three-momentum 
vanishes, we have 


P= (p, £1), P2 > (—p, £2) 


and hence 
_ tpl(E + Es) 
. Fy, k, 
puting this value of ua in 
An? BE E> 


Im f(a > a) = — (—UaFa/ 167°) 


E 
finally we get 
k; 

I > a) = —d4 

m f(a —> a) in? 
Where k is the total magnitude of momentum of incoming particles in C.O.M frame. 
This form of the unitarity condition is known as the optical theorm. There is a pretty 
consequence of the optical theorem that tells much about the pattern of scattering at 
high energy. The scattering amplitude f may be expected to be a smooth function of 
angle, so there must be some solid angle AQ within which |f|? has nearly the same 


value (say, within a factor 2 ) as in the forward direction. The total cross-section is 
then bounded by 


1 1 
aa > f [fd > SIf(a + 0)PAM > 5|Im fla > a)PAD 
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Using the optical theorm then yields an upper bound on AQ 
AS 32n?/k’ ox 


As we shall see in the next section, total cross-sections are usually expected to ap- 
proach constants or grow slowly at high energy, so above equation shows that the solid 
angle around the forward direction within which the differential cross-section is roughly 
constant shrinks at least as fast as 1/k? for k > oo. This increasingly narrow peak in 


the forward direction at high energies is known as a diffraction peak. 


B. Boltzman H-theorem The same argument that led from the unitarity condition 
StS = 1 to our result 


Im Mag =—7 f d85* (pp — Po) Meal 
also allows us to use the other unitarity relation S.St = 1 to derive another result 


Im Mag =—7 f d85* (pp — Pa) [Maal 


Putting these two equations then yields the reciprocity relation 


[485 pa — pa) Maal? = f 864 (pp pa) Mol? 


or in other words 


(gees = (j0038 


This result can be used in deriving some of the most important results of kinetic theory. 
If P,dx is the probability of finding the system in a volume da of the space of multi- 
particle states ®,, then the rate of decrease in P, due to transitions to all other states 
is P, [ d8dI'(a — 8)/d, while the rate of increase of P, due to transitions from all 
other states is [ d3PsdI'(8 + a)/da; the rate of change of P, is then 


dP, dV(8 a) _ dV(a — 8) 
aan porno Pa f dp iB 


It follows immediately that | P,da is time-independent. (Just interchange the labelling 


of the integration variables in the integral in the second term in above Equation, On 
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the other hand, the rate of change of the entropy — [ daP, In P, is 


d 
- a daPein Py =~ f do f d6(nP+1) 
: [p, te a) pte >| 





Interchanging the labelling of the integration variables in the second term, this may be 


d a Ps \ aV(8 > a) 
-5 [ daPoin Py = fda f dsPsim (2) are . 


Now, for any positive probabilities P, and Pg, the function Pg ln (P3/P,,) Partial- Wave 


written 





Expansions satisfies the inequality 


Qa 


P 
P3In (2) > Ps — Py 


The rate of change of the entropy is thus bounded by 


d 


_ dl'(8 > a) 
dt 


daP, In P, > : da / dB [Ps — Pal 7 
(ey 


or interchanging variables of integration in the second term 


i di(g—4a) dI(a- 8) 
-5 [dorsinr.> [da fasrs| da ~—~— dp 





But using the unitarity relation 


jgtee (ghee 


tells us that the integral over a on the right-hand side vanishes, so we may conclude 


that the entropy always increases: 
d 
—— | daP,InP, >0 
dt 


This is the famous Boltzman’s H-theorem 
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3.1 Path Integral Formulation 


3.1.1 Introduction to path integral formalism in Quantum 


Mechanics 


3.1.1.1 Action Principle and Sum over all paths 


Quantum Physics and QFT can be formulated in three ways: 
1)Matrix Mechanics(Heisenberg 1925) 
2)Wave Mechanics(Schrédinger 1926) 
3)Path Integrals(Feynman 1942-1948) 


The first two versions, whose equivalence was demonstrated very early by Schrodinger, 
Dirac, and others, are covered in every quantum mechanics course or textbook; the 
third is the focus of this thesis. It was created by Feynman and is based on Dirac’s 
work . For a long time, this formulation was thought to be excessively complex and 
worthless for practical reasons, thus it was left out of textbooks. The Feynman path 
integral, on the other hand, is easier to comprehend conceptually and needs a con- 
siderably less drastic departure from classical physics concepts (to which we are all 
accustomed) than the conventional quantum mechanics with its operators and state 
vectors in a Hilbert space. When it comes to the practical side of the various formu- 
lations, it is undeniable that Schrodinger’s wave mechanics is the quickest and easiest 
approach to get a solution for many problems: for example, when computing the sta- 
tionary states of a particle in a given potential. However, in many-particle situations, 
the Schrodinger equation is considerably less important, and Heisenberg’s matrix me- 
chanics is much more important (e.g. when diagonalization the Hamiltonian matrix 
in a certain subspace). Similarly, in field theory, the path integrals technique has be- 
come an essential tool: The Feynman rules for non-abelian gauge theories were initially 
obtained in this way, and numerical attempts to address the strong interaction field 
theory are based directly on the (Euclidean) path integral representation of Quantum 
Chromodynamics. Furthermore, this approach differs from other formulations in that 


it may be used to a wide range of physics problems. This ” unification” of disparate 
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subjects not only aids comprehension, but also suggests new techniques that have been 
successful in other domains. One caveat: Except for certain circumstances (Euclidean 
formulation, i.e. for imaginary times Brownian motion Wiener integral), the Feynman 
path integral is not rigorously defined mathematically. Numerous attempts and for- 
mulations have been made to offer a mathematically solid foundation . We won’t be 
concerned about this in our thesis because, in most situations, mathematical rigour 
lags well behind physical intuition. Quantum mechanics was established by Feynman 
in a manner similar to classical mechanics. As a result, let’s start with a quick rundown 


of the most important parts of classical and quantum mechanics. 
Classical Mechanics 


A system (A point particle in one spatial dimension) with coordinate(S) q(t) is de- 


scribed by the Lagrangian 
L(q(t), a(t), t) 


It’s dynamic development, that is, its trajectory (or “path”) from an initial point 
q(ta) = da to an end point q(t») = q@ proceeds in such a way that among all possible 


paths the one is selected that makes the action 


ty 
s= fata. 4.0 (3.1) 
ta 
minimal (or more precisely:external): 
6S =0 


From that equation of motion (Euler-Lagrangian equations) follow 


dOL OL 
—— — — 3.2 
dt 0q oq ° 2) 
and define the hamiltonian as 
(p(t), a(2),t) = [P-4—L4.4.lama (3.3) 
Hamilton’s equations 
oS OH 
q Op 
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can also be obtained from an extremum of the action 
to 
15 =6 f atlp-4— H(p(0).4(0).t) = 0 (3.4) 
ta 


by independent variation of q and p. 


Quantum Mechanics In classical mechanics Nature seems to “compare” different 
paths and then to select the path of extremal action. This is also the case in optics, 
as long as the wavelength of light is much smaller than the typical dimensions of the 
system. However, if both sizes are comparable, we observe the typical interference 


behavior of waves, e.g. diffraction by a double slit as depicted in Fig 3.1. 


The different ways now play a crucial role, since the scattered wave is composed of two 
parts 


® x= cikli 4 eikle 


where k is the wave number . The intensity at the detector is proportional to |®|? and 


shows the typical diffraction minima and maxima. 


yr 


double slit screen 


Figure 3.1: Diffraction of waves through a double slit. 


It is now a crucial experimental fact that even matter particles (like electrons or neu- 
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trons) show diffraction effects when passing through slits or scatter from crystals. This 
happens when their de Broglie wavelength is comparable to the dimensions of the ob- 
ject to be imaged (for electrons, this was demonstrated for the first time by Davisson 
and Germer 1927). 


The quantum mechanical description of this effect postulates a probability amplitude 
which is a superposition of the two possible amplitudes (for the passage through the 
respective slit): 

@= 6,4 ® 


and the probablity of finding the electron at the detector is , 
W ~ |o/ 


in other words one has to sum coherently over the various alternative ways an event 


can happen. 


Let us imagine that we now drill more and more holes in the screen — until it no 
longer exists — and that we put more and more screens in the space between source 


and detector which we also treat in the same way, then we have 


(a,b) = Ss" ®; 


all path from a to b 


Of course, it remains to specify more precisely how this sum over “all paths from a to b 
” has to be performed and with which weight each path contributes. The answer to the 
first question leads to functional integrals, which we will treat in the next chapter while 
the problem of the weighting is solved in a heuristic way: Since Planck’s constant h 
has the dimension of an action and determines the quantum effects it is not surprising 
that the rule is 
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t 


Figure 3.2: Quantum mechanical propagation of a particle from a to b. The 
classical path is denoted by a red line 


All paths contribute with the same absolute magnitude but with different phase: The 
phase of every path x(t) is its classical action S[x(t)| divided by h 


®(a,b) = > const.e(Ol/" (aco) 


all path from a to b 
With this prescription the classical limit (formally: h — 0) can be performed imme- 
diately: If all actions S>> h this leads to immense oscillations of the exponential 
function, that is, to each path with a positive contribution, there is an adjacent path 


which gives a negative contribution and therefore deletes it. There remains only the 


special path for which this doesn’t happen, viz. whose phase is stationary 


Oca = 0 


and only one trajactory , contributes to the sum of all paths. 
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3.1.2 Path integral representation of Evolution Operator /Dynamical 


propagator 


In the following all quantum operators are denoted by an ”hat”. we also assume that 
the Hamiltonian 


A A 


H=T+V,V =V(a) 
is time-independent and consider matrix elements of the time-evolution operator 


Ui (ite. bar iast = (2 10 (to, ta) ie — (2p eel a) (3.6) 





As we all see, the time-evolution operator contains all information about the system 


that we need. 
Feynmann’s Path Integral 


The path integral representation of U is based on the decomposition of the time interval 
t» — tg into N single, infinitesimal small intervals 


th — tp 
N 


so that 
N 


el -ilte—ta) H/h] = jim eich /h (3.7) 
k=1 


This is also called Trotter’s product formula. Now it is fact that 
el-ie(P+V)/h] _ p[-teT'/A] . el-teV/R]+O(€?) (3.8) 
i.e. for infinitesimal times the non-commutativity of the the quantum mechanical 


operators can (almost) be neglected. If we use this expression in the Trotter’s formulae 


and insert the identity operator 


+00 


(oe) 


between each factor, it is only a small step to obtain the path integral: Since for a local 


64 


3.1. Path Integral Formulation 





potential V(z)|x) = V(a)|x) is valid, we obtain 


+00 


Ute eas he) = slim dx ,dxy...drn_1 (2 jeer tw-1) e eV (en-1)/h 
oo ee 
; (en-1 ect! ty-2) ciV (ena) (3.9) 





: (2, ian wa) ei (wa) /h 


we now need the matrix elements of e(-*!/") between position eigenstates 
2 
~ieT /h es ok a a (-+#) ip;-(xj41—25)/h 
Davie 2 ar _ D;€ € (3.10) 





Here we have used that T' = f?/(2m) is diagonal in momentum space. The integral in 


(3.10) is an extended Gaussian integral which can be derived from the basic integral 


a 2 7 
G(a) = i) aie = /% Rea > 0 
i a 


by completion of the square: 


+00 42 
/ da.e(-a0? +be) _ zel#) ,a#0,Rea>0 (3.11) 


(oe) 


with this result we get 


ae) = Cais: 312 
73 Qrich . ( ) 





and if, for notational simplification, we set 
To = La, UN = Xp 


we have 


; m N/2 -++oo 
Ui, tet, be) — sim (—) / dx,dxq...dtn_1 


N-1 (3.13) 
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The product of exponential term is the exponent of their sum and we obtain 


N/2 ea 
=a) is dx 1 dx Susie dxn—4 
(3.14) 


ier |m (x gone 
1-2; 
“exp ey (3 ( | — V (z;) 





UO (Litt a) im ( 
= 








The last line is a symbolic shorthand for the limit 
N>w,e7>0 sEeN =ty —ta 


of the (N-1)-dimensional integral. It is indeed the action which appears in the argument 
of the exponential function in this limit because the discrete sums approximate the 


Riemann integral 








(3.15) 


Note that the short -time approximation 
exp[—ie(T + V)/h] = exp[-ieV /h] - exp[-ieT'/h] + O (€) (3.16) 


leads to an action in which V(x; + 1) appears instead of V(a,;). In a similar way the 


more symmetric form 


exp[—te(T + V)/fi] = exp[-ieV /(2h)] - exp[—teT /h] - exp[—ieV /(2h)] + O (€?) (3.17) 


would lead to a slightly different discrete action. In the simple case under discussion 
(one particle in a scalar potential) this does not lead to any consequence: The differ- 
ences are of order € and vanish in the continuum limit « + 0. The situation is different 
if one considers velocity-dependent interactions because the velocity is represented by 


(x;41 — x;) /e. More on that under the heading ” Ordering Problem”. 
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The expression (3.14) gives an exact prescription how to sum over "all paths from a 
to b” : One discretizes the time and approximates each possible path from a to b by a 
polygon path. Integration over an intermediate point x; corresponds to the summation 
over all paths from x;_; to x;41, because with fixed x;_1, 7j;+1; a new path is generated 
by a different value of 2;. There is no integration over the endpoints because they 


are fixed by the boundary conditions. An extraordinary aspect of the path integrals 





tase § bg eae Me 


Figure 3.3: Discretization and summation over the quantum mechanical 
paths. The red path arises from a different value of 7; 


is that they deal with ordinary numbers and not with non-commuting operators! This 
is an enormous advantage compared to the usual formulation of quantum mechanics. 
The price one has to pay for that is hidden in the complicated definition of the path 
integral: It is an infinite-dimensional integral, a so-called functional integral. It can be 


shown that the usual Schrodinger theory follows from Equation 3.14 . 
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Phase Space Path Integral We obtain another form of the path integral if the p 


-integration of the Gaussian integral is not performed: 


+00 


im dx\dx- cate dxn—4 / 


a | 
N-oo gs none 


3 N 2 
LE Le Vie D; 
0 ED [Ds J = 2 -veva)} 


j=l 


a / “= Dix(t)Dp(t) 
7 a(ta)=ta Qrh 


ti dp, dp2 dpn 
Qrh2rh ~° 2rh 





U (Xp, tp; Xa, ta) 





exp] 5 faint) - 20) — Hote) atop} 
(3.18) 


Some comments about the “phase space” or Hamiltonian formulation: 


1. The argument of the exponential is again the classical action divided by the action 


quantum, but this time formulated in coordinates and canonical momenta . 


2. The “integration measure” dxjdp;/(27h) is different from the Lagrangian formula- 
tion; for distinction we therefore use the notation D’x(t). This is descriptive and easily 
remembered as dividing the phase space into cells of the size h = 27h (as in statistical 


mechanics). 


3. There is an additional integration over dpy and no boundary conditions for the 


momenta. 


4. The interpretation is more difficult: Is this a sum over trajectories in phase space? 


Classically, however, a single point in phase space determines the trajectory. 


5. While in the Lagrangian formulation the paths are continous and only the veloc- 
ities are discontinuous (see Fig. 3.3), in the phase-space path integral the paths are 
discontinuous as well. Therefore caution is needed when performing canonical trans- 


formations and other manipulations . 


6. Klauder has shown that by introducing a regularizing factor 


; 1 ty -2 
Jim DxDp exp {Flea} - exp -=, | dt (2a? + =) 


the phase space path integral can be defined unambigously. 
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The Ordering Problem If a classical Hamiltonian is given, in which products of 


coordinates and momenta occur, the question arises how to quantize such forms 


”No rule based on the correspondence with Classical Mechanics can resolve such am- 
biguities, since the latter arise from the non-commutability of operators, which in turn 
is tied to the finite and non-zero character of h... In all cases of practical interest, 
one must conform to the following prescriptions: ... Having put the Hamiltonian func- 
tion into the form >>, p; fi (q1... dn) one replaces... (this) term... by the symmetrized 
expression 7 i aif at i ibi 

How do these ambiguities show up in the path integral? Obviously not in the symbolic 
notation as used in Eq. (1.30) and (1.46)! Actually, they appear in the innocent- 
looking problem how the Hamilton function is treated in the discretized form of the 


path integral, e.g. 
a) outer averaging: }[H (p,;,2;) + H (p;, 7;-1)| 


71) 


b) inner averaging: H (p;, 5 


c) right hand: H (p,, x;) 
d) left hand: H (p;, x;-1). 


If there are no products of x and p in the Hamiltonian, then - as discussed before — the 
various procedures give identical results and one can chose the simplest formulation. 


In the general case one can show that for the treatment of products of the form 
A= pe” (3.19) 


in the Hamilton function, there exists the following correspondence between the above 


path integral prescriptions and the operator ordering: 
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Since #,p are hermitean operators, only the symmetric forms a) and b) also give her- 
mitean operators. Rule b) is called Weyl’s quantization rule or ” mid-point rule” and 


is obtained automatically if Eq. (3.19) is considered as Wigner transform 


+00 


A= Aw(z,p) = i 


—co 


dy (x- Ale + 3) elpylh (3.20) 


of the corresponding quantum mechanical operator A. The reverse transform is 


- Oe ata’ . : 
A ‘) = Bay ip (w—a')/h 21 
(2 le i mh w( 2 ye eal) 


If we apply this to each factor in Trotter’s product formula, we obtain 


ha FOO iy. os . : ; 
al x;) -| =r Cae (3 “2,)) ets (ej41—3)/h 
ee, 








Co 


=1—ieHw ((0j41+0;)/2,p;)/ft-... (3.22) 
+00 dn. ; ; 
Pj v€ jy t Xj Ay Ge ATR 
~ rad ——H ‘ DpALj+1— Xj 
i ma | hi w( 2 w)]e 


which leads to a generalization of Eq. (1.45) for arbitrary Hamilton functions and 
which contains the mid-point rule. The ordering prescription becomes important for 
curvi-linear coordinates or for quantization in curved spaces. A less exotic application 
is the case of a charged particle in a magnetic field because the minimal substitution 


for the electromagnetic field leads to the velocity-dependent Lorentz force. 


3.1.3. Perturbation Theory and Retrieval of Interaction Pic- 


ture Time Evolution Operator 


We now consider the (for simplicity again: one-dimensional) motion of a particle in a 


general potential V(x), i.e. Lagrange functions of the form 


In general the path integral cannot be evaluated anymore,A widely used method for 


weak interactions (potentials) is perturbation theory, which we may obtain from the 
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phase space formula (3.18) by expansion in powers of the potential: 





x(t')=a! D'rD i U p 
EM Gree =| L +f 4 ——— 
Creat) byte, oan exp | = o\ pi — 5 


| 
Me 
ai 
S| =. 
See 
GS 

= 

® 

a 

8&8 

cee 


(3.23) 


The first term is the free propagator which is given by 


+00 d é 9 
Up (Xp, ty; Gide) = 7 a exp E (> ‘ (Xp a a5) _ al 





a (3.24) 
= Le (Zp — La) 
oninT > | on Me Fe 
Put 1 =2,% =xandt —t=T 
+00 d : 2 
he gh: = ap ‘a r _ ee —_ 
Uole. oa) = ie: aah xP . |» (x’ — x) a (t 0] \ (3.25) 


Let us now consider the next term of the perturbation expansion 3.23 


x(t')=a! D'rD i i p t! 
oat, = P ey pes. 
UE e eh) = osc. aR xP if do (0 F / drtV (x(T)) 


In the discretized form the last factor is 








and thus 
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= +00 +00 
dx, dpr, dpn 
: ([ 2th cae ann’ a) 


Similarly as in the evaluation of the free propagator we write the exponent as 


eed iN oy 
a Pp 
i b» (pi — Pi41) Li + PNT’ — pit — € d z 


and perform all x; -integrations (except x;). This gives (NV —2) momentum 6 -functions 


balk: a oe dp N-2 
OG te) = sim Oy Pe TI ea (27h) 


—oo 


-d (pi ae 6 (po — p3)...6 i — ps) 6 (pj41 — Pj+2)-..6 (pw—1 — PN) 


N 54 
a Pi 
“exp {i (pj — P41) 2j + pnz’ — pia] —€ 2 HI 
These enforce p; = po = ...pj-1 = pj and p41 = pjt2 = ...PN-1 = pn, So that we 


obtain 





rat ae _ oes a ! 
Ui @,8) = sim Oy ey (25) (Onh)? exP 4 [(pj — pw) 23 + pwx — p;z] ep - 


(3.27) 
ee [ie (Bi , Px w=) (3.28) 
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If we now write € = x;,p = p,; and p' = py, the result is 


+00 co dr! A 12. 
iat a) = = Jim of dEV (E f- ser YF w-o9'- Faw -a|} 
i 


=Uo (x ,¢/ 3€,t+Je) 


+00 : 2 
[How tg |€-2- Eiie . 
a, 


=Uo (€,t+je;x,t) 





(3.29) 


It is now possible to perform the limit which gives a Riemannian integral over the 


intermediate time 7, (7; = t + je) 


t! +00 
U; (a’,t/;2,t) = / dr / deUp (2! ,ts,7) VQUd(E, 32,8) (3.30) 


Read from right to left this contribution to the perturbative expansion can be depicted 


as a Feynman diagram as shown below 


time time 





space space 


Figure 3.4: (a) Feynman diagram in 1st order perturbation theory for the 
time-evolution operator in Eq.(3.30).The solid blue lines represent the free 
propagator between two space-time points, the dashed red line the interac- 
tion of the potential V . Intermediate points are to be integrated over. (b) 
Feynman diagram for the 2nd order, see Eq. (3.32) 
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Similarly one proceeds in higher orders 


1 pee =e' py «Dp ‘78 P t! 
Z tat, =—_— _ i 
UG Gra, 6) i i ah XP |e / do (0s F / drV (ax(rT)) 


By induction one can prove that 


(f" seven) = fie i. or a drtjV (a (11) V (a (72))...V (a (74)) 


which - with a similar calculation - leads to 





+! +oo 
Use ttt) = / ar f COE Gr) VIGO Ne Tat ae ce. B81) 
t —oo 


The graphical representation for 7 = 2 


t! +00 
Us et at) -| anf dEQU  (x", t'; €2, T2) V (E2) 
—00 (3.32) 


T2 +00 
i, an | dE, Up (€9, 72; 1,71) V (1) Uo (Gir e52) 


is shown in Fig. 3.4b . Insertion of the recursion relation (3.31) into Eq. (3.23) gives 


an integral equation for the time-evolution operator 
i t! +oo 
Ula tat) Une tiat)— 5 far f aedo le’ ts6) VOUETA2) B33) 
t —cO 


This is the usual equation for the time-evolution operator in the interaction picture 


as one can see if one writes above Equation in the operator form 


. t! 
G (t',t) = Op (t,t) — 7 / dry (t', 7) VU (r,t) (3.34) 
t 


Defining 
U (t', t) =: Uo (t', t) Uz (t’, t) 
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one obtains 


3 t! 
Up (t’, t) Uz (t’,t) = Uo Ce i dr Up (t',T) VU(r,t) (3.35) 
t Ne 


=U (t! ,t)Uo(t,7) =Uo (t’ t) Ug (7,t) 


Here the composition law for the free propagator Up has been used. Multiplying from 
the left with Uj! (t’,t) one obtains 


Ur (t’,t) =1- ai dr Ut (r, t)VUo(7, t) Or(7,t), 
; i ee Eee 
=VY (t/t) 
which is identical with the equation of motion for the time-evolution operators in the 
interaction picture. From the derivation of Eq. (3.33) it is clear that this equation is 
also valid for time-dependent potentials. However, if the potential is time-independent 
then the integral term on the right side of Eq.(3.34) is a convolution integral, since 
the evolution operators then only depend on the time difference. According to well 
known theorems for the Fourier transformation we then immediately obtain for the 
time-independent Green function 
+00 


G(E) = i dtG(te”/", with G(t) = —Lo(itt 0) (3.36) 


(oe) 


which simply is the operator identity 


i] 1 I “ 1 
V 


E-H+i0+ E-A)+i0t+ E-M+i+t E-A+i0+ 





The prescription that the singularity in these operators has to be regularized by adding 
+7x an infinitesimal positive quantity in the denominator follows from the definition 
(3.36) in which the energy £ must have a small, positive imaginary part in order that 


the time integral converges. 


75 


3.1. Path Integral Formulation 





3.1.4 Greens Function as Path Integrals and Time ordering 


in Path Integrals 


Up to now we have discussed the path-integral representation of the time-evolution 
operators U in different formulations. Now we want to consider 2-point Green functions 


of the form 


Gap (tata) = (0]T [An (tr) Bu (t)]|0) (3.37) 
where | 0) is the exact ground state of the system, 
Aner Ane 
is a quantum mechanical operator in the Heisenberg picture (as is B,(t)) and 
- An (1) Bu (t2)| = © (ti — tz) Au (t1) Bu (ta) + © (t2 — tr) Bu (t2) Au (th) 


the time-ordered product of these two operators. Such objects play a central role 
in many-body physics and in field theory, as it is straight-forward to evaluate them 
in perturbation theory and all important quantities (ground-state energy, expectation 
values of operators, S matrix elements) can be calculated from them. If we insert into 


the definition (3.37) a complete set of position eigenstates 


t= f ealay(a = fave ™ayalen™”. 


then we obtain 
cilt'/h ‘) (¢ ett! /ha- An (t1) Bu (t2)| eft a) 


F (4 eee 0) 


The individual factors in this expression may be rewritten in the following way: The 


(0 ‘) = eiBot'/h (Q Il”) = eiBot' [Pays (7) 


can be expressed by the wave function of the ground state and the ground-state energy 





Gap (ti, te) = [ eadd (0 








(3.38) 


first one 
eilit! /h 








- similarly, also the last one. Therefore we now concentrate on a path-integral represen- 
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tation for the middle factor and first assume that A = B = ¢ and that t, > ty. After 


inserting two complete sets of position eigenstates we then obtain 
(a e iit! [hae (t,) 2 fa (ty) en a) = ( ! a) 


= | eadoanae (4 1| iA! “IM a) Can lem ty—te "| as) (ae Jew to— afl @) 


Now we may use the path-integral representation of the individual time-evolution op- 





aq Lie 5 i es ‘K os 
ei t'—t1)/fip pil (t-te) /hg tft (tot) 














erators and obtain in the Hamilton formulation 


(d ett he (t1) 4 (is) eittih @) 


a(t')=q' D'«Dp i t! 
= | dqd dripz — H 
/ 71 woe fo Oonh exp a T [px (x, p)] 
a(ti)=q1 D'xDp 17 ty 
exp / irlpi — H(2,9))} 
—_ Qh h ta 


POD EDD { Ef? 
. exp i; dr|pz — H (x, p) \ 
ee. 2Th h t | 


Now we can concentrate the exponents into one exponent and all integrations can be 














combined into one functional integral if one integrates additionally over the endpoints 


41,42. For t; > ty the result therefore is 


(qe "TF [ar (ta) bx (to)] "| q) 
a(t!)=q' D'xDp i v : (3.39) 
=, oa x (t1) x (te) exp =f dt|pa — H(x, p)| 


)=4 





For tz > ¢,, a similar calculation gives the same expression, i.e. the time-ordering 
operator is automatically built into the path-integral representation! This is 
not quite so surprising as it first seems, since the time ordering is necessary because the 
Heisenberg operator do not commute at different times — the path integral, however, 


deals with ordinary numbers. 


Eq. (3.39) is easily generalized to any operator A = A(#,p),B = B(&,p) and then it 
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becomes 


(a eile (ne An (t:) Bu (t2)| aun a) 


x(t!)=4q' Dl 
7 i ee (p (t1) ,# (t1)) B (p (ta) ,« (t2)) exp {pst co] 





(3.40) 





Let us return to the formula (3.38) for the Green function in which we now can insert 
the above result. However, the appearance of the ground-state wave function w9(q) is 
unsatisfactory since the exact ground state is unknown in general. We would like to 
have an expression available in which the ground state is generated at the same time. 
This is possible by ” filtering” out the lowest lying state of the system. Indeed, in an 


expression of the form 
/ 
(q 


only the ground state contribution survives if one formally approach by doing a Wick 





emo 4) = yee (q) 7 tEmt!/h mm] O|n) eBnt/y* (q) 


rotation and get it in the form of a Eucledian Path Integral 
t—ioo, t!' > —ioo 


and get 


@ 





euniogs™ a) re ase (q’) Wa (qe Bol+lD/* (91 OQ) 


fas . 3.41 
= lim (d tea q) (0/00). oo 





This is due to the fact that for n > 0, E,, > Ep, so that the contributions of the excited 
states decay more rapidly than those from the ground state. If we use above Eq. (3.41) 


for the Green function of the operators A, B we obtain 





santtuty = tg (CET [Ae Bet oo 
AB (1, %2) = Um 











Fi (3.42) 
Pipe (a! eH —8/n a) 
or as path integral form 
__ f D'rDpA (x (t1) , p (t1)) B (2 (te) , p (ta)) el" 
G gp (ti, te) = oe i DiaDpeSteallt (3.43) 
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Note that in the expression (3.43) all normalization factors cancel! The unphysical 
limits t > ioo, t’ + —ioo are quite natural for Euclidean Green functions where the 
time t is replaced by —i7 everywhere. Instead of oscillations one has damping every- 
where and the path integral ist well-defined mathematically. Therefore one frequently 
investigates Euclidean Green functions and transforms back - if possible - to real time 
by an analytic continuation. It is clear that in order to project out the ground state 
by the boundary conditions t + ioo, t’ — —ioo one doesn’t necessarily have to go 
to infinity along the imaginary axis but that it is sufficient to do that along a suitable 


ray in the complex time plane 


The operators A, B can be expressed in terms of the fundamental operators #,p, say 
by a Taylor expansion. However, it is not necessary that in the matrix element the 
products of operators <,p all are taken at the same time: One can define n-point 


correlation functions in quantum mechanics as 


G (ty te .t;, tj4a os by) = (0 \T [tu (t) ay fH (t;) Du (t541) So PH (tn)]| 0) (3.44) 


which takes the form 


PD aDpe Gee aa) ap (te) eiSlx,p]/h 


G (ty bya tj, jay ee oe ty) if D'xcDpeiSlzl/h 


(3.45) 





in the Path Integral Formulation. 
The full set of n -point functions can be obtained from a generating functional 
i 


Z| J, A= [ PeDvexp {pst + - ieao + K(pp(0)} (3.46) 


by functional differentiation with respect to artificially introduced sources which are 
set to zero thereafter 

of am Z|J, K] 
Od (tr) 0. dd (€;) 0K (t541) ... 0K (t,) 20, 0] 








J=k=0 
(3.47) 
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In general the generating functional cannot be evaluated exactly. However, if one can 
split up the Hamiltonian 

H = Hy) +V(@,p) 
in such a way that the generating functional Zo|J, AK] is known then one can give a 


closed expression for the full generating functional: First one also splits up the exponent 


in the path integral as 


/ dt [pi — Ho(w,p) + J(t)x(t) + K(#)p(t)] - : atV (2(t), p(t)) 


and then one uses the fact that x(t) can be generated by functional differentiation w.r.t. 
the source J(t), p(t) by differentiation w.r.t. A(t). With this procedure the exponent 
of the ’perturbation” V can be taken out from the path integral and we obtain the 
representation 


a 


Z|, K] = i D'sDpexp | . / WV (2(#. 040) exp | 50+ 5 / ays bel + K(Ho(0)} 


=p -j fav Ge "TKD) | fos Fe 





(3.48) 


Of course, this is only a formal solution and in most cases one has to expand the 
exponential function in powers of the perturbation in order to perform the functional 
differentiations. This produces the perturbation series for the generating functional and 
therefore for alle n -point functions. We will discuss this in detain in the next sections 
in Field Theory Path integrals. 
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3.2 Field theoretic treatment of Path Integrals 


3.2.1 Functional Calculus 


3.2.1.1 What are Functionals? 


Mathematically Given a manifold W representing (continuous/smooth) functions f 


(with certain boundary conditions etc.), and a functional F defined as a map 


F: 4-98 or F: 4 ->C, (3.49) 


It is evident that in dealing with dynamical systems, we are dealing with functions of 
continuous variables. In fact, most of the times, we are really dealing with functions of 
functions which are otherwise known as functionals. If we are considering the motion 


of a particle in a potential in one dimension, then the Lagrangian is given by 
; 1. 
Dae pe gine” — V(a) 


where x(t) and «(t) denote the coordinate and the velocity of the particle and the 


simplest functional we can think of is the action functional defined as 


Sigh [ 2 arid 


Note that unlike a function whose value depends on a particular point in the coordinate 
space, the value of the action depends on the entire trajectory along which the inte- 
gration is carried out. For different paths connecting the initial and the final points, 
the value of the action functional will be different. Thus, a functional has the generic 


form 


FIfl = if deF(f(2)) 


where, for example, we may have 


Sometimes, one loosely also says that F(f(x)) is a functional. 
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3.2.1.2 Functional derivatives 


The functional Derivative or the Gateaux derivative of a functional F (linear functional) 


= gods, a 
_=fs 5f (a) (x) (3.50) 


is defined as 





Eg “rif + ev] 





where v is an arbitrary function and ev is defined as variation of f Equivalently, from 
the working point of view, this simply corresponds to defining 


bF(f(2)) P(f(@) + 0(@ — y)) — FUF(@) 


—_—~ = lim 


f(y) «0 € 





It now follows from above equation that 


Of (zx) 
éf(y) 


The functional derivative satisfies all the properties of a derivative, namely, it is linear 





= (rz — y) 


and associative, 











6 _ Fils] , Sls 
afta) I+ BUD = Sra) + afta) a 
5 _ OF If] dF >If] . 


It also satisfies the chain rule of differentiation. Furthermore, we now see that given a 


functional Ff], we can Taylor expand it in the form 


Fif] = [ acPy(a + [ew dx2P; (41, £2) f (x2) 


+ [ow dz dx3P, (x1, v2, £3) a (x2) a (x3) Se 


where 
Pot) = FUF(@))| pe)=0 
éF (f (x1)) 
Of (x2) f(x)=0 
Lae (f (@1)) 
2! Of (v2) of (xs) 


Py (histo) 








P, (1, Lo, L3) = 





f(x)=0 
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and so on. 


Functional differential: The differential (or variation or first variation) of the func- 
tional Ff] is 


OF | f;v] = | Geomaas 


Heuristically, v is the change in f, so we formally’ have v = df, and then this is similar 
in form to the total differential of a function F' (fi, fo,..-, fn), 





Or 
dF= 5 OF df; 
i=1 


where fi, fo,..-, fn are independent variables. Comparing the last two equations, the 
functional derivative 5F'/d f(x) has a role similar to that of the partial derivative OF'/0;, 


where the variable of integration «x is like a continuous version of the summation index 


2. 


3.2.1.3 Functional Integration 


Whereas standard Riemann integration sums a function f(x) over a continuous range 
of values of x, functional integration sums a functional Gf], which can be thought 
of as a “function of a function” over a continuous range (or space) of functions f. 
Most functional integrals cannot be evaluated exactly but must be evaluated using 


perturbation methods. The formal definition of a functional integral is 


fowler ={_-f_enlle@) 


However, in most cases the functions f(a) can be written in terms of an infinite series 


of orthogonal functions such as f(x) = f,H,(x), and then the definition becomes 


fovtras ff eUut. oT 


which is slightly more understandable. The integral is shown to be a functional integral 
with a capital D. Sometimes it is written in square brackets: [Df] or D[f], to indicate 


that f is a function. 
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Examples: Most functional integrals are actually infinite, but then the limit of the 
quotient of two related functional integrals can still be finite. The functional integrals 
that can be evaluated exactly usually start with the following Gaussian integral: 


fel-3 (y)dedy+ J J(e)-F(@)de/ 7 f) = eth f I(x) K-M(ay) I(y)dedy 
rar Seg) 





By functionally differentiating this with respect to J(x) and then setting to o this be- 
comes an exponential multiplied by a polynomial in f. For example, setting K(x, y) = 
6(a — y), we find: 





[F@OI@ETIOO® Dis ee ae 
f eff f@)F@)aatDFl = K"(a,b) = la — bj? 


where a,b and x are 4-dimensional vectors. This comes from the formula for the 
propagation of a photon in quantum electrodynamics. Another useful integral is the 


functional delta function: 


peires z)dei 1) f) = d[g = [1s (3.52) 


which is useful to specify constraints. Functional integrals can also be done over 
Grassmann-valued functions ~(x), where w(x)w(y) = —w(y)w(x), which is useful in 
quantum electrodynamics for calculations involving fermions. It is worth noting 
that Functional Integrals and Path integrals are same thing and only dif- 
fers in the terminology used by Mathematicians(Functional Integrals) and 
Physicists(Path integrals). 
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3.2.2 Path Integral formulation of Correlation Functions(Field 
theory) 


3.2.2.1 PI formulation of the two-point correlation function 


A further step toward a more convenient formulation of QFT is obtained by getting 
rid of the boundary conditions. What we are really interested in, is the set of time- 
ordered correlation functions among fields at different space-time points, evaluated on 
the ground state of the system (also know as Green’s functions). Indeed, we will see 
how , starting from these correlation functions the LSZ reduction formulae allow us to 
extract the probability amplitudes for physical scattering processes, i.e. the elements 


of the S matrix.The simplest of such correlations is the two-point function 


(2|7 {n (e1) dt (v2) $2) (3.53) 


where |Q) denotes the ground state (of the interacting theory) and bx the operators 


associated to the fields in the Heisenberg’s picture. 


In order to find a PI formulation of (3.53), consider the following functional integral, 


with generic boundary conditions at t = +7 


hig = | Do(e)6 (21) 6 (aa) FH (3.54) 


The main functional integral [ D¢(x) can be broken up as 


[Pee = [Paola f ; 


0 
1 
0 
o(29, 


namely introducing two specific boundary conditions at times x? and x9 (in addition to 
those at the endpoints —T and T ), but integrating separately over the corresponding 
fixed-time field configurations ¢;(#) and 2(Z).Consider first the case T > 23 > x? > 
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—T. Under this time ordering, the exponential term in (3.54) can be decomposed as 


0 
ci ltp tt el(o) _ iene aLl(?) x lad d4xL(4) e cite tall) 


After the decomposition of the main integral, the extra factors ¢(21) and (x2) in 
(3.54) become ¢; (% 1) and @2 (#2), and can be taken outside the inner integral. We 
thus obtain 


id= / D1 (Z)Dbo(#)or (21) bo (#s) 














x (or eo 2) ; (2 ei (#8- x ey ) ; (2 fe —iH( («2 +7) ) 61). (G3 (3.55) 
SS 
if], £()dx4 09 ed 
[Doe 2 Ms pea ee (pee rte 


O(T) = or and ¢(—T) = ¢; .The main integral then factors into three pieces, each 
being a simple transition amplitude as shown above. ,Note that the times x? and 23 
automatically fall in order. We can turn the field ¢; (%,) into a Schrédinger operator 
using $s (£1) |¢1) = 1 (1) |¢1). The completeness relation [ D¢1 |¢1) (¢1| = 1 then 
allows us to eliminate the intermediate state |¢,). Similar manipulations work for ¢2, 


yielding the expression 


0 


(or je) bs (x2) e tH (22-#1) 4 (x1) ew tH(steT) 1) ; 
Most of the exponential factors combine with the Schrodinger operators to make Heisen- 
berg operators. In the case x? > x$, the order of x; and x2 would simply be inter- 


changed. Thus expression 3.54 after similar steps becomes 


(bp |e" T {by (21) bu (2) $e" | b7) - 


This expression is almost equal to the two-point correlation function. To make it 
more nearly equal, we take the limit T — oo(1 — ie). Just as in Section 3.1.4, this 
trick projects out the vacuum state |Q) from |¢,) and |¢y) (provided that these states 


have some overlap with |Q), which we assume). For example, decomposing |¢,) into 
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eigenstates |n) of H, we have 
oF? Iga) = > (e7 |n) (0 | ba) 7-s00(1-se) —? (2 | Ga) e OPE), |) 


As in Section 3.1.4 , we obtain some awkward phase and overlap factors. But these 
factors cancel if we divide by the same quantity as (3.54) but without the two extra 


fields @ (a1) and ¢ (x2). Thus we obtain the simple formula 


a, DH (1) 0 (20) exp [i [Zn dae] 
Wm. 
T-00(1—ie) | Dé exp i ee del] 





(Q [Pou (#1) bn (#2)| Q) = (3.56) 


This is our desired formula for the two-point correlation function in terms of functional 
integrals. For higher correlation functions, just insert additional factors of @ on both 
sides.Once we do that In general we get 

_ f Doe S49 (a1) +++ 9 (an) 


Ci (Bis Benes Dy) =O () Ol Gs) > Oy) = J Does (3.57) 





For time-dependent correlation functions, the time-ordering operator T is included. 


This is known as the n -point correlation function also called n point greens 
function in Canonical QFT.Further are gonna use it to derive our S -Matrix in Path 
Integral formulationin the last chapter.It should be noted that Correlation functions 
are also called simply correlators which can be interpreted physically as the amplitude 


for propagation of a particle or excitation between two space time points say y and x. 


3.2.3 What are Generating Functionals in QFT? 


The generating functional Z[J] is a key object in quantum field theory - as we shall 
see it reduces in ordinary quantum mechanics to a limiting form of the 1-particle 
propagator G'(z,2’ | J) when the particle is under the influence of an external field 
J(t). However, before beginning, it is useful to look at another closely related object, 


viz., the generating functional Z(.J) in ordinary probability theory. 
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Recall that for a simple random variable ¢, we can assign a probability distribution 


P(@), so that the expectation value of any variable A(¢) that depends on ¢ is given by 
(4) = f doP(o)AC) 
where we have assumed the normalization condition 
[P(o)=1 
Now let us consider the generating functional Z(J) defined by 
ZisS / dbP(¢)e® (3.58) 


From this definition, it immediately follows that the” n -th moment” of the probability 
distribution P(@) is 








O”"Z(J) 
= (o") = | dbP(d)d" = 3.59 
m= (6) = f aPC =O) (3.59) 
and that we can expand Z(J) as 
= cea 
Z(J) = 3+ — | dbP(o)o" 
x m | (3.60) 
1 


so that Z(J) acts as a “generator” for the infinite sequence of moments g, of the 
probability distribution P(¢) 


For future reference it is also useful to recall how one may also expand the logarithm 
of Z(J); we write, 


(3.61) 
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Now suppose we make a power series expansion of W(J), i.e., 


wan=>- =CnJ” (3.62) 
n=0 ~ 


where the C,, known as ”cumulants”, are just 


arw (J) 


Cn = Bn 








J=0 


The relationship between the cumulants C;, and moments g,, is easily found. We notice 


that a = 7 
OW(J) OmMZ(J) 1 0Z(J) 











OD. SOF = fh Ou: : z 
WJ) 8 (1dZ(J)\_ [1 Zs) 1 AZ(J) 
OF OINZ OF 7° (VA OF Fe OF 


and so on. Continuing in this way, and using eq. of C,,, we can find the {C;} in terms 
of {g;}, by expanding In Z[J] = In {exp ae = ond” }, and comparing coefficients of 
J”. We then find that 


=%9 

C2 = 92 — gi 

C3 = 93 — 39192 + 29%, 

C4 = 94 — 395 — 49193 + 129297 — 6gt 


and so on. Finally, note that in simple probability theory, we can make a ” Legendre” 


transformation to change the dependent variables; if we write 


W(J) =1(¢) + Jo 


then in terms of the new function I'(¢), we have 


and we can then rewrite everything in terms of '(#) and its derivatives. We will 
see in what follows that all of these functions of simple probability theory can be 


put into correspondence with functionals in QFT. As we will also see, it is useful 
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and illuminating to explore the analogies between probability theory, field theory, and 


thermodynamics/ statistics mechanics. 


3.2.3.1 Generating Functionals in Field theory 


Suppose we are now dealing with a quantum field ¢(x) in D dimensions. Then the 
analogue of the generating function in eq.(3.58) for ordinary probability distribution is 


a generating functional Z[J], given by 


Z| J] = N f DoAle: J = Ojetf Prd(@) oa) 
(3.63) 
= N f Debits s?7ebve 


and we see that the analogue of the probability distribution P(@) is the ” partition 
function” Z[J = 0, ¢], ie., the normalization factor N in the above equation is given 
by 

Nt= j= | Doei/*5141 = / DoA(¢ 


A(¢) = = ei/hS[¢) — ei/h f dP aL ($l #)) 


(3.64) 


so that the analogue of the probability distribution P(@) for a random variable to take 
the value ¢ is now the amplitude A(¢) = e’/"*!@l for a field to take a given configuration 
(x). We shall see later what we must write to deal with probabilities (as opposed to 
amplitudes) in QFT, using a kind of ”double path integral” technique (recall that in 
ordinary QM one discusses probabilities using the density matrix, which also can be 


understood in terms of a double path integral). 


Note also that a convergence factor has been left out of eqs.(3.61) and (3.62). Just 
as in ordinary QM, we need to specify which configurations are being integrated over, 
and with which boundary conditions. In ordinary QM, this is done by rotating the end 
points of the time integration to +ioo, or by adding a term to the exponent to force 


the Green function [ Dge'/"*\4 to pick out the ground state hence we write 


ZL =N f DpA(d; I= oe #a[seieorree 
(3.65) 
_y i: Dei! S 2 2[L(0(@))+Ja)o(e)-+ie0?(@)] 
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I do not go through the demonstration that this gives the desired result, since it is 
similar to (but more complicated than) the demonstration for ordinary Q.M.. In any 


case, the key point here is that 
Z|J] = (0 | 0).1(2) 


i.e., Z[J] is the vacuum expectation value for the theory, under the influence of an 


external source J(x). 


3.2.4 Time ordered Correlation functions from Generating 


functionals 


The time-ordered correlators/correlation functions in QFT is defined by putting Time 


ordering in equation 3.57 as 


: e/FSIA 6 (21) ---b (tn 
Gn (21,°-* ,2n) = (0|P {6(m) 6 (n)}] 0) — Le ee ats 





(3.66) 


where the central term is the actual defination and the last integral term is the Path In- 
tegral form of the correlation function the denominator is added to normalize the result 
(we cannot in general assume that the function is normalised); then the corresponding 


power series expansion is then 
Zi = Son few Z »f dene, ECON ACH EEC 
n=O nr 


and the G;, is related to Z[J] by 


2 [J] 


Gn (X4,° ee it) = (—ih)" GF (x1) -6J Ca) 








J=0 


Noting that Z(J) in ordinary probability theory is just the expectation value of the 


function e?7, we see that in scalar field theory, we have 


z= (0 lf jonieee 0) 
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which when expanded out gives 


ZI] = (010) +> f a?x(0l0(x)I0) 10a) 


+a faa, ip dary (0|F {6 (a1) 6 (a2)}] 0) J () J (02) + 





3.2.5 Path Integral Green’s Functional of a scalar Quantum 
Field 


The transition from many-body quantum physics to field theory is a transition from 
a finite number of degrees of freedom to an infinite number. We will perform this 
transition formally thus omitting in a first step all subtleties like renormalization, even 
ignoring the question of the sheer existence of an interacting theory.In the beginning we 
consider for simplicity a system of particles of mass m which is described by a neutral 
scalar field ¢(x). Let the Lagrange density (Lagrangian) be 


L = = (J, 80"® — m6’) — V(®) (3.67) 


1 
2 
where V(®)denotes the self-interaction of the field. If one requires the theory to be 
renormalizable (in 1 time and 3 space dimensions) this interaction can only be a poly- 


nomial of the field up to degree four; a typical example is 


V() = A6*/4! (3.68) 


which is of great importance in superconductivity (Ginzburg-Landau theory) and par- 
ticle physics (Higgs mechanism). As in the quantum-mechanical case we consider the 


transition matrix element 


(®, le *#Ur-#) 





®;) (3.69) 


of the time-evolution operator between field configurations ®; and ® and H = f[ d®?xH 
is the Hamilton operator of the systems. The classical Hamilton density follows by the 


usual Legendre transformation as 
; i ee 3 bbe 
H=nr®-L= 57 + 5 (V2) + 5m &* + V(®) (3.70) 
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where 1 = OL/0® is the canonically conjugated momentum (density) of the field. We 
divide the (finite) space volume L? into N small cells with volume v = L*/N and the 
time interval t,t; in M intervals of length o. If we use the completeness for the fields 
Phi(t,x) = ®j where, 1=1...M 1,j=1.. .N at each time and space point Z 


/ de! |!) (| = 1 (3.71) 


then we obtain 


(8, Je“#r“4)| @, 








_ M-1 1 M | .-ievH| @M-1 
)= st. f TL do! (o le @! ) 


.(O} Je] wo) b 


ore we have set ®// = ®;(x,) and ®) = ©, (x,) as abbreviation. Exactly as in the 


—ievH 








quantum mechanical case finds for small times 


141 | —ievH 
ay" fe 





ia [ Prexp [iv' (oot? — ®') — ievH'] 


l 


If we write Pi = U7}, we obtain 





(®, le *#ler-*) 





» (et -x)\ (3.72) 
2/2On2 pat) haf i at f d's [rb — (x, 0))} 


This is the expression for Path integral of a sclar field in Hamilton form . By 
completing the square (no — 7?/2= —(n — 6)?/24 $/2) we can transform the 7 
-integral into a Gaussian integral and we then obtain the Path integral in Lagrange 


form as 


(8, je*(er-4) 





;) = const. » f DoG e's (3.73) 
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Note: Again we have the (classical) action in the exponent of 3.71 because 
[ at fi +6 +(V)? 1 ee? V(o) 
‘ me o2 5" 
tf 1 1 tf 
-| at f ds 52,00" — se = via) ay at f d'xe (®(x),0,®(x)) = S/O] 
tt tj 


However, this representation of the transition matrix element between field configura- 
tions is only an intermediate step as this is not the relevant quantity in field theory. 
Actually, all observables in field theory can be obtained from the n-point 


functions or Green functions given by 3.64 
Cae (0|7 [@(e1)...& @,)} 0) 


(as can be seen below in the example of scattering of two scalar particles). Here |0) 
is the exact ground state of the theory (the trivial ”vacuum”), 7 the time-ordering 
operator and ®(x) are the exact field operators in the Heisenberg picture. What we 


need therefore is a path-integral represention for the Green functions. 


Fortunately, this already has been in the previous sections for case of quantum me- 
chanics and is not hard to generalise it to Quantum Field theory . The ground state 


can be projected out by the unphysical limit t; > 100,t’ — —ioo and one has 


“1, He Jf D®® (x1)... ® (an) exp(iS[®}]) 
Go@aesetn) = “oly J Dd exp(is[@) 


tf 4— coco 





(3.74) 


The set of all n-point functions of a scalar field can be obtained from the functional 


differentiation of Generating functional Z|] of the scalar field which is 
Z[J] = / D® exp (‘sia +4 / teJ()®(2) (3.75) 


Hence G,, can be expressed as 


ol 5” 4\J] 
CE (=) 6J (x1)...6J (an) Z[0] 





(3.76) 





J=0 


As one can see the normalization of the functional integral cancels. Eq. (3.73) can also 


be considered as functional Fourier transform of exp(iS[®]) and thus the artificially 


94 


3.2. Field theoretic treatment of Path Integrals 





introduced source J(a) can be seen as variable conjugate to the field ®(2). In most 


cases the generating functional can be calculated only for vanishing self-interaction V : 


: ZolJ] = / D® exp f it d‘x (5a.0ar — sme? - 0)] 
= / DO exp E if d(x) (0? — m2) ®(x) +i / ode) %(a)) 


In the second line an integration by parts has been performed under the assumption 
that boundary terms vanish. The functional integration in the above Eq. can be done 


since the exponent 


: i d‘xd*y®(x) Ko(x, y)®(y) +i / d'x I(x) P(x) 


s 5 [aay (®Ko + J) (x) Ky (x,y) (Ko® + J) (y) — 5 [ aratys(e) Kye.) I) 


is quadratic in the fields. The Gaussian integral over ® thus gives 


a 


Ae -5 / dndty I(x) Ky Ua, nau) = 70] exp | (J, Kz") 


(3.77) 


a 
2 
where the constant does not depend on the source. In addition, we have used a conve- 
nient short-hand notation. The inverse of the kernel 
Ko(,y) = (-0) — m*) o°(a —y) 


is calculated most easily in momentum space as - due to translational invariance —Ko 


only depends on the difference x — y : 


A(e,y) = Ky '(e,y) = ‘| SaA(Hete 


The equation 


then just becomes an algebraic equation 
(2 — m2) A(k) = 1 
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with the solution 


Apr(k) = k2 — m2 + i0+ (3.78) 





This is the Feynman propagator in momentum space. Of course, one has to specify 
how to treat the pole at k? = m? this can be achieved by a more careful treatment of 
the limit t; + ioo,t, + —ioo in the path-integral representation (3.74) for the Green 


functions, or simpler by introducing a convergence factor 


into the path integral. By combining this term with the mass term in the Lagrangian 
this leads to the replacement 


ne SS oe a0” 


(i0* is a short-hand for a small, postive imaginary part te,¢ > 0 which at the end of 
the calculation is set to zero.). From Eq. (3.75) one thus finds that the 2-point function 


(or the propagator) in lowest order is given by 


GO (x1, 22) = iAp (21, 22) (3.79) 


The generating functional with interaction can be obtained formally from Z| J] by 
functional differentiation. As functional differentiation w.r.t. the source J(x) generates 
the field ®(a) we thus have 


Z| J] = [> exp -i fatev Ga) exp ji [ay (Lo + J(y)®(y)) 


ue 1 / dav Gl Zol J 


(3.80) 
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3.2.5.1 Machinery of Generating Functionals:A Functional Equation for 
the Generating Functional 


The representation of the generating functional as a functional integral offers many 


possibilities for transformations known from usual integral calculus. For instance, we 


may perform an infinitesimal shift of the integration variable 


O(x) = v(x) + ef (x) 


As the Jacobian of transformation is obviously is “1”, thus we have 


ARIE [re {1 tie fof exp (‘ste +i f atzs(ayo(a)) }. 


Since the function f(a) is totally arbitrary, the (functional) integral over a total (func- 





tional) derivative must vanish (as in ordinary integral calculus): 


0= [Pesaow (‘ste +i [ atzs(a)o(a)) 


= i[ De { ae + 112) exp (iste +i [ ates(a)o(a)) a 





dp(a 


We may pull out the terms, which contain y in the curly bracket, from the functional 
integral if we replace y(x) by 6/(idJ(x)) (acting on the generating functional) in these 


terms. Then we obtain 





{Bey sal : a) 20 =0 


which for the scalar theory explicitly reads 


1G +m”) ID + Vy! (ar) — a) 20 =0 (3.82) 


This so-called Schwinger equation therefore connects different n -point functions 
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3.2.6 Greens function as Feynman Propagator 


In this section, you will make your first acquaintance with path integral. In quantum 
theory, path integral is a method of constructing the propagator. So we begin with the 
propagator. Recall the Green’s Function ,mathematically Green’s function G(2; x) of 


an operator L is defined as 
LG(a;2 ) = 6(a;2) 


this method is very useful to solve equations like 


since if we know the Green’s function G(x; 2) of L, the solution is simply 


(zr) = f Gas!) ae’ 


One can easily verify it notice 


= [tee f J (e da = f 5(0:0") Ia) = JI(a) 


Where the second equality is due to the fact that L consists only operation on x and 


is totally irrelevant to x’. 


In Quantum Mechanics , we have the Schrodinger’s Equation 


O 1 
s— + Vv? =() 
f Ot = 2m | Iv) 
All within the bracket can be be viewed as an operator, and thus possesses a Green’s 
function a F 
feet __w2 K eghial\pa jawed ace : 
i+ sv | Ct Oe Oe) (3.83) 
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with this propagator , we can calculate the state of anytime after from an initial state 


|w(t’)) = [Ke x;t', x')|w(t))dadt (3.84) 


therefore , the initial state |q)(t)) propagetes to the final state |w(t’)) with the help of 
K(x,t;x’;t’). As a consequence,K(x,t;x’,t’) is called the propagator. 


Feynman suggest the propagator takes the form 


K= / Dre's (3.85) 


3.2.7 Perturbation theory in Path Integrals 


The relation 


Vii em 1 i‘ dav Ga) ZolJ (3.86) 


between the Generating Functional with and without Interaction is purely a formal one 
and comes to use only when we expand the Interaction in terms of coupling constant 


of the theory. For this consider the expansion of form 
Z[J] = Zo[J] {1 + Awy[J] + wo[J] +...} 


where 


i SAE jee (aia) ml 


Expanding the differentials gives 


wih -+ i. ds {3(—1)’Ap(a,2)Ap(a, 2) 


+6(—i)3 / Oy d'yAr (0,11) J (gh) Me(@,2)Ar (0,92) F (2) (3,87) 





4 
+ i) fant yat pty, [[ 4. (x, yj) J (yj) 


j=l 
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The different terms can be ordered according to the number of sources J and represented 
graphically: Each propagator Ap(x,y) is symbolized by a line leading from y to x and 


each interaction vertex by a point into which 4 lines are going in or out as shown below 


Y1 Y3 


(A1) (B1) (C1) 


Figure 3.5: Graphical representation of first-order perturbation theory w 
for the generating functional in ®4-theory 


In second-order perturbation theory one finds 
iz : 
wal J] = swild] + Hel] 


where the first contribution on the r-h.s. contains disconnected graphs in which two 
parts are not connected by a common line (see, e.g. Fig.3.6 below where the first-order 


graphs (A1) and (B1) are taken together without connection). 


Figure 3.6: A disconnected graph in 2nd order perturbation theory for the 
6*-theory 
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3.2.8 Connected and Amputated Greens Functions and their 


Generators 


We have seen how the generating functional or partition function or vacuum expec- 
tation value Z[J] corresponds to the generating function Z(J) for simple probability 


distributions. Let us now continue this analogy, and write 


Z[J] = ek Wl (3.88) 


We now proceed , by defining what we call ”connected” Green functions (the name 


will become clear) as 


1 wl] 


GO (21,...,2n) = (ey a 5J (#1) +++ 5F (an) | p29 


(3.89) 





and 


Wy] = 5 (5) 5 far. .aege Goce Tsay. 2650) 


n=1 


Composing eqs. (3.89) and (3.90), we see that 


Coie ea = Z[ Go (Bij es) 


3.2.9 Vertex Functional of Scalar Field Theory 


So far we have expanded two different functionals as infinite series, i.e., the vacuum en- 
ergy generating functional Z[.J] (the partition function), and the generating functional 
W|J| for connected diagrams (analogous to the free energy in statistics mechanics), 
related by 

WJ] = -thln Z[J] (3.91) 


with expansions 
aU (5) aL f aan 21-529) J @) (3.92) 
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and 


Now we introduce a 3rd expansion, in terms of a new generating functional I'|¢]; this 
however, is not a functional of the external currents, but of the field ¢(x) itself. We 


write 


T[¢] =o fan fatsr, (iid gin NO (a) 22.20 (Eas) (3.94) 


so that 
é’T [J] 


59 (21) +» 06 (an) 





| BP es gerne ere (3.95) 





o=0 


Now of course the way to change variables between the J(z;) and the ¢(z;) is by a 


Legendre transformation, with form given here by 


rg] = Wi] - y de I(x) $(2) 


Thus now we now have the pair of relations 











OW | J] as 

0F() ly (3.96) 
TL _ sey) 

dg(x) J 





Vertex functions like T, (v1,...,%,), turn out to be central to the approach of QFT 
to the dynamics of real physical systems - they are used in everything from scattering 


theory to particle physics to transport theory in condensed matter systems. 

3.2.10 Relation between Connected Greens function and Ver- 
tex Functionals 

Here we will simply establish the relationship between the connected Green functions 


go (%1,..-,2%,) and the proper vertices [,(x%1,...,%p). Let’s first do this for the 


simplest case, that of the 2 nd-order correlators go (41,2) and Ty) (41, x2). From their 
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definitions, and from eq. (42), we have 





























ow | J] do (71) 
GO (x1, 22) = — 1? ae 
2 (1,2) bJ (#1) 6F (x2) | 59 6J (x2) | 55 
TJ] OJ (x1) 
[5 (21,22) = == 
aia ae XC (CT) Pema Ce) PO 
from which we immediately find that 
c do (x1) 
dyG (x1, Pegs) = =Wd(aj=—2 
i YY (71,4) 2(y 2) 50 (x3) (74 2) 
whose Fourier transform is 
—f? 
gi) _— = 
2 (PP) T'3(p, —p) 


i.e. the two functions are inverses of each other. 
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4.1 Scattering Matrix (S-Matrix) construction in 


Path Integral formalism 


Now we have developed all the necessary tools in order to formally construct our S- 
Matrix in the Path Integral as well analyse its Unitarity. We start by giving a rough 


sketch of how we are going to proceed. 


4.1.1 Sketch of Theoretical construction 


Before we step into the path integral form of S-Matrix , we give an overall sketch as 
well as some rules of our construction here. Beware, that the rules are not necessarily 
accepted or respected by all the academics.The construction will start from analyzing 


the general path integral of scalar field from some stable state to itself (usually it is 
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vacuum state |Q)) 


(Q | Q) =n [ Dye {—5i f ¢ [a — m?| o} =a (4.1) 


where JN is a normalization factor. The last equality comes from the property stable, 
which means it will not change under the situation of free field, which will determine 
the normalization factor. Then, we construct a quantity called generating functional 


by adding an arbitrary potential J 


Step 1 
WolJ] = (2) 2), =n [ Dpew{-i f 5¢ [a? — m?] e—Jo} 
=exp{—5i [ tars} 


The above generating functional includes merely the scalar field without any other 


(4.2) 


actual interaction terms (the J term is an auxiliary term with no physical significance). 


If the Lagrangian includes the interaction terms £;, , the generating functional will be 


wis] =exp{-i f em 22) wats 


cool fel pon 2 0] 


As is indicated, this expression deals with some stable state. Now we hope to gener- 


Step 2 


(4.3) 


alized it into other ordinary states, like from state a to state G. This will give us the 


well-renowned S-matrix (scattering matrix) in P.I formulation 
Step 3 


.) 


$= (8 0) = STlyaa = em {i fu (8*— me) 5 hv] (4.4) 





J=0 
where Yas is the asymptotic field corresponding to the initial and final state. The last 


equality 4.4 is nothing but the famous LSZ reduction formula(but now in Path Integral 


formulation). 
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After arriving at the S-matrix, the theoretical construction will come to a conclusion. 
The S-matrix, which is directly observable in the scattering experiment, may be used to 
calculate the scattering cross-section. As previously indicated, we will not incorporate 
any components from canonical quantization into our development of the path integral 
formalism. As a result, we are unable to utilise notions such as state, which are 
frequently used in quantum theory. In this scenario, we’ll need to employ some field 
setting as a state (like asymptotic field). However, this implies that there must be 
a subtle identical connection between field and state, which canonical quantum field 
theory rejects. In actuality, the difficulty caused by the lack of a field operator and 
consecutive state only emerges when creating an S-matrix, thus stage 1 and stage 2 
are unaffected in any case. The path integral formalism, as a complement to canonical 
quantization in our philosophy, should be free of notions in canonical form. As a result, 
we'll do our best to avoid using terms like field operator and matching state, or at the 


very least push them as far down as feasible. 


4.1.2 Free Scalar Field theory 


This section will introduce the free scalar field theory, i.e. stage 1 mentioned in the 
sketch. Our main goal here is to derive the generating functional of the free scalar 
field. First we need a vacuum state as the initial and final state — just like the start 
and end point in path integral of Quantum Mechanics. However, the vacuum is not 
so straight-forward in Quantum Field Theory. We hereby define the vacuum state 
|Q)as the configuration of field with the lowest possible energy, which is the general 
starting point of the term vaccum. Now, we express the probability of the free field 
from vacuum state to vacuum state in terms of path integral. The probability should 


be unity since there are no interactions. Therefore 


W = (2/9) =¥ [ Doc | -3i f (o-m) o} = (4.5) 


where we now know the normalization factor NV. We already know the Lagrangian. 
The scalar field is real here, and it won’t be miserable to generalize this into complex 
situation. However, This formulation does not tell us anything except for a normaliza- 
tion factor, since it is trivial - from vacuum to vacuum with no interaction. Now, we 


want to add some interactions. To achieve this, we introduce an auxiliary potential J 
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and the Lagrangian becomes 


L=-p (0 as m?) ptJlyp (4.6) 


and the path integral becomes 


WU] = 2] 9)I, =N [ Dees | -5i fv (— m2) p- 249} (4.7) 


Now, we perform the path integral. As what we have done in Quantum Mechanics, 


expand the path integral according to the classical path. The classical path yo satisfies 
(0 _ m?) Yo = J, ie. Yo = [acs 
where A, is the Green’s function satisfying 
(0? — m?) Ap (a; 2’) = —6 (2; 2’) 


Warning: minus the Please about hand s side the be the sign above equation! very 
careful in right h of tl Consequently, we perform translation y + y+ %p. In path 
integral, a definite path yo acts just like a constant (expand the path integral according 


to the definition and you will see why). A brief calculation will give 


WJ] =exp{—5i [ tars} N | Doew {-3i [ o(e-m?) o} 
=exp{—5i [ tars} 


4.1.3 Scalar field theory with Interaction 


(4.8) 


This section introduces the self-interaction theory of scalar field. As usual, the in- 
teraction is implemented as an additional interaction term Lin, in Lagrangian. The 
interaction discussed here is assumed to satisfy the following condition It is the self- 
interaction, i.e. Ling = Lint [y] It has the form of polynomials. Similarly, the construc- 


tion starts from the generating functional. The generating functional of an interaction 
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W[J] = [ Peextis) = [deem {i f c+ emlel} 
= [ Peexr {i f enlel} exw {i fc} 


Next, we expand the first exponential (the one with interaction) according to definition 


wi = [Pes i fen lh] ex uifes 


Now, recall our first identity 





theory is 


(4.9) 


L 0 


Fa Wolll = | Devle) exptis} 


Therefore, each v(x) in Ling [py] can be replaced with 6/idJ(x), as long as Lint [y] is a 


pure polynomial of y. In this case, we can write 


vin= for’ fou ah Fe 
-£3[+ fo EAD] frvenlef a 


where expression Ln; [6/idJ] represents the the expression where all y in Ling is re- 
placed with 6/idJ(x). The terms after the replacement no longer contain y and thus 
can be moved out of the path integral. And we find now that the path integral gives 


the generating functional Wo|./] of free field. Transform back into the exponential form 


Woe {-i / Lim ba \ Wold] (4.10) 


which is the result of second stage of our construction. The interaction of other kinds 


and we get 


of fields (except for vector field in gauge theories) can be performed through similar 


Manner. 
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4.1.4 Scattering Matrix 


Scattering matrix, or S-matrix, is one of the major and core concept in both Particle 
and High Energy Physics - since this is the most straight-forward way to verify the 


form of interaction. 


As we have already seen that the S-matrix is defined as an infinite dimensional matrix 
with each matrix element S53. being the probability of the state evolving from a to £, 
Le. 


Spa = (8 | a) 


Note: As a matter of fact, the S-matrix is not a typical ” matrix”. It is just a quantity 
with each one of them can be labelled by two indices. Also, in canonical quantization, 
S-matrix is actually an operator which, if it is finite dimensional, can be expressed 
as a matrix (but in Quantum Field Theory it is not). Therefore, when you see some 
concept which is a matrix, it might not mean to be a typical matrix, but some operator 
or quantity of two indices. In Quantum Field Theory (as well as Particle Physics), we 
are especially interested in the following situation Two or more beams of free particle 
are generated and targeted at some point, interact with each other at the limited 
region near the point, then leave the interaction region, become free particles again 
and finally detected by the particle detectors. The key point in the above picture 
is that the interaction only happens at a restricted region, and the particles are free 
outside the region. This approximation can be pretty accurate if the region is selected 
wisely so that the particles are separated by a considerably long distance. And it is 


prevalently implemented in many particle colliders. 


The above picture suggests us write the field in the path integral as the addition of an 
"interact field” y and a free asymptotic field yas , i.e. perform the following variable 


translation in the path integral 
PIF Pt Pas 


where the asymptotic field y,, satisfies 


(0 — m?) G0 
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Therefore, the S-matrix with auxiliary field J supplied is 


Sold] = [Pee \—i fe | — m?*] e— Je ~ Jae} = exp {iPasT} WolJ] (4.11) 


Similar formula also holds for situation with interaction. Then, just as what we have 


done in the last section, expand the exponential as Taylor series 


oar n) 
Recall we have our second identity 
(2? — m?) = 7 Wold] = Ja) Wold 
i dJ (zx) 


This inspires us that we can replace every J in the Taylor series with (0? — m?) 6/id J. 


Hence 
1 I, Bc odie tales! fe ; gf asa liso 
S\= dS iPas (O = ees WJ] = exp § ipas (0 SH) WIJ] 
n=0 


To obtain the S-matrix, simply set the auxiliary field J as zero 


5 = exp {ives (0% me) 25 bw 





(4.12) 


This is the so-called LSZ reduction formula, which gives the S-matrix with inbound 


and outbound state designated by asymptotic field y,, from the generating functional. 


4.1.4.1 S-Matrix in Terms of Correlation Function 


Recall that the Correlation function that we studies in the last chapter can be expressed 
here as 


1 1 6 
seh r5 Ger 7 id J (21) oF OD ay! J=0 





which is the correlation function, or n -point (Green’s) function. The origin of this 
name comes from the similarity of its formulation in canonical quantization and hence 


this name may look strange here. It is not so straight-forward about the n -point 
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function, and thus we hope to reduce the formula. To do this, let us first evaluate the 


2-point function 








1 6 1 6 1 
T (1,82) Sepia | if rsh 2% 


stay | O (rai) Jes {Fi f Jed} 


Geant ( i aed) Wald 





J=0 


= iAp (x1; £2) 
J=0 





It is nothing other than i times the ordinary Green’s function. Next, 3 -point function 


at (fe) 


= 0 
J=0 


T (1, %2,%3) ~ Ap f ApJ 








And 4-point function 


T (1, 22+ X3, 24) = iA (21; £2) iA (23; La) 


+ iAp (21; 23) iAp (223 24 





) 
ait iAr (2; £3) iAp (v1; 4) 
=(1A VB. 4+ N3)2084+4 200 3)0~4) 


where we have used (1 ~ 2) to represent 7 (21,272). Therefore, we can see the regu- 
larity — the n -point function is either 0 or some combination of two point functions. 
Specifically, we have by induction The odd-point function is 0. The even-point function 
is 

F (Giese) = SS (pi ~ pz)... (Pen—1 ~ Pon) 
which means the summation over all the permutations p. This is the so-called Wick’s 
Theorem - transforming the n -point function into the combination of 2 -point function. 


It can be easily shown that the S-matrix can be written using correlation function as 


$= >> SJ] oe) TL  -) + (er, (4.13) 


where x takes over 71,...,%,. This formula is the so-called perturbative expansion 


of the S-matrix in Path integral formulation. 
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4.1.5 First glance at Infinities 


This section will show you how the infinities emerge in Quantum Field Theory through 


a resonated toy of physicists —y* theory. The interaction term of y* theory is 


1 
Le = age" (4.14) 


Recall that the generating functional of interaction theory is 


WJ] = exp {oi fe Fa] \ Wold 


Suppose the interaction is weak enough for us to neglect all terms of order higher than 


one. Therefore, the generating functional up to the first order is 


wf] = (1 f on Fal }) Wold] = ( 2 ral Bl} Wold] 


Notice that in this case the four functional derivatives are at the same point x. There- 


fore, a calculation similar to the previous section will show 


to. 
i dJ(x) i dJ (x) 





WolJ =a) SA ee ire 


We see that the right hand side contains the value of Green’s function at zero, which is 
clearly infinity. With higher order terms which contains even more derivatives, there 
might be even more infinities. To deal with this infinities, we use a mathematical 


technique called renormalization. Looking back to where the infinities come and we 


find 
— (= i Ar(eiy)J(u}av) 


makes no sense, since the functional derivative requires that A(x; y) be a function, 
but it is essentially a distribution (or generalized function). But things are not so bad. 
The Green’s function is a distribution which usually has a function correspondence, 


and for most of the physical equations that we care, it has. 


Nonetheless, the function corresponding to A(x; y) has no definition on x = y. There- 


fore, the above expression hopes to get something undefined. This makes things much 
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easier. If it is undefined, the most straight-forward way is to place an additional defi- 


nition. 


The renormalization procedure does exactly this stuff. Usually, we will impose some 
conditions (such as define the physical mass as the square-root of minus zero-point of 


two-point convex function) so that the value A(x; x) can be defined. 
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5.1 Problem of Quantum Gravity 


Reconciliation of quantum mechanics with gravity may well be the most profound theo- 
retical problem left unsolved by twentieth-century physics. Both general relativity and 
quantum mechanics became well established in the first quarter of the last century, yet 
their clash continues to this day. One might simply declare the problem too difficult, 
and move on to more modest goals. But, there are glimmers of hope, both from new 
approaches to the problem, and perhaps more importantly, through a better under- 
standing of the nature of the underlying difficulty. There are a number of sub-problems 


to that of quantizing gravity. A list includes: 


1. The problem of UV divergences and nonrenormalizability: what structure cures 
an apparent infinite proliferation of perturbative ultraviolet divergences, and needed 


counterterms? 
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2. The problem of singularities: how can a quantum theory resolve singularities of the 


classical theory, such as are found in black holes and in many cosmologies? 


3. The problem of observables/time: how can one formulate gauge invariant observ- 
ables in the quantum theory, given that diffeomorphisms are gauge transformations in 


the classical theory? 


A. The problem of unitarity: what dynamics unitarizes amplitudes, in the high-energy 


regime? 


5. The conundrums of cosmology: this includes the problem of defining infrared finite 
observables in inflationary cosmology — an extreme example is the measure problem; 


or, perhaps we need an alternative to inflation. 


Historically, there has been a lot of focus on problems one and two, and they have 
strongly motivated approaches to quantum gravity such as string theory and loop 
quantization. However, there are increasing indications that problems three through 
five are both more profound, and are linked; they seemingly drive more deeply at the 
heart of the central problem of quantum gravity. Indeed, there has been a growing 
appreciation of the fact that there are deep long distance issues in quantum gravity, 
which do not appear easily resolved by short-distance modifications of the theory. This 
suggests that a successful theory must incorporate more radical departures from local 
quantum field theory (QFT), relevant over even very long distance scales, in certain 


circumstances. 


5.1.1 Approaches to Quantum Gravity 


The quest for quantum gravity can be separated into three main lines of research. The 
relative weight of these lines has changed, there have been important intersections and 
connections between the three, and there has been research that does not fit into any 
of the three lines. Nevertheless, the three lines have maintained a distinct individ- 
uality across 70 years of research. They are often denoted “covariant,” “canonical,” 
and “sum over-histories,” even if these names can be misleading and are often used 
interchangeably. They cannot be characterized by a precise definition, but within each 
line there is a certain methodological unity, and a certain consistency in the logic of 


the development of the research. Despite the fact that general relativity has been 
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present for almost more than 100 years,and have passed the subsequent test of time as 
well still most physicists have indeed agreed that it cannot be the ultimate theory of 
gravity.This is because of astounding success of the Standard Model and the Quantum 
Field theories that lie underneath the standard model.These theories which includes 
both abelian and non abelian groups has been verified experimentally to a level of 
precision that the Generel relativity cannot match. Another straightforward reasons 
include the breakdown of GTR in the vicinity of black holes and the singularities given 
by the theory. It is also a fully classical theory, whereas all other physics theories 
are quantum mechanical in nature.On the observational side, we now have extremely 
precise verifications of general-relativistic effects in the solar system, as well as rather 
compelling evidence for strong-field predictions like black holes and the ”big bang.” 
While some unproven conjectures remain, such as cosmic censorship, the development 
of new mathematical techniques has given us a fairly complete qualitative understand- 
ing of the theory, while computer advances have enabled us, at least in theory, to make 
quantitative predictions to any desired order of accuracy. The final argument, however, 
remains valid: despite a great deal of effort (and some results), we still lack a suitable 
quantum theory of gravity whose classical limit is general relativity. This is most likely 


today’s most significant unresolved topic in theoretical physics. 


Historically There are three main ways of quantizing the classical gravitational field 


theory. 


e The covariant line of research is the attempt to build the theory as a quantum 
field theory of the Perturbative fluctuations of the metric over a flat Minkowski 
space, or some other background metric space. The program was started by 
Rosenfeld, Fierz and Pauli in the thirties. The Feynman rules of GR were la- 
boriously found by DeWitt, Feynman and Faddeev in the sixties. t’Hooft and 
Veltman, Deser and Van Nieuwenhuizen, and others, found increasing evidence 
of nonrenormalizability at the beginning of the seventies. Then, a search for an 
extension of GR giving a renormalizable or finite perturbation expansion started. 
Through high-derivative theory and supergravity, the search converged success- 


fully to string theory in the late eighties. 


e The second method is the canonical approach in which one introduces a family 
of spacelike surfaces and constructs a Hamiltonian ( eg ADM formalism of Gen- 


eral Relativity ).In this formalism the Hilbert space carries a representation of the 
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operators corresponding to the full metric, or some functions of the metric, with- 
out any background metric to be fixed. The program was initiated by Bergmann 
and Dirac in the fifties. Unraveling the canonical structure of GR turned out to 
be laborious. Dirac, Bergmann and his group, and Peres completed the task in 
the fifties. Their cumbersome formalism was drastically simplified by the intro- 
duction of new variables: first by Arnowit Deser and Misner in the sixties and 
then by Ashtekar in the eighties. The formal equations of the quantum theory 
were written down by Wheeler and DeWitt in the middle sixties, but turned out 
to be too ill defined. A well defined version of the same equations was successfully 
found only in the late eighties, with the formulation of LQG. Although many peo- 
ple favor this, the division into space and time seems to be contrary to the whole 
spirit of relativity. Also it is not clear that the concept of a spacelike surface 
has any meaning in quantum gravity since one would expect that there would 
be large quantum fluctuations of the metric on small length scales. Further, we 
shall want to consider topologies of the spacetime manifold that do not permit 


any well-behaved families of surfaces let alone spacelike ones. 


e Third one is the “Path Integral Approach” also known as the sum-over- 
histories line of research is the attempt to use some version of Feynman’s func- 
tional integral quantization to define the theory. The idea was introduced by 
Misner in the fifties, following a suggestion by Wheeler, and developed by Hawk- 
ing in the form of euclidean quantum gravity in the seventies. Most of the 
discrete (lattice-like, posets,. . . ) approaches and the spinfoam formalism, 
introduced more recently, belong to this line as well For these reasons we prefer 
the path-integral approach though it too has problems concerning the measure 
and the very meaning of the integral.But still this method seems to be the most 


fundamental one as we have seen its success in the other quantum theories. 


Now we will briefly discuss all three above mentioned major approaches and we will 
construct and discuss the behaviour of the S-Matrix in all these formulations and 
analyse the Unitarity and analiticity ,thus giving us an idea of which one of these 
approaches seems most promising, since Unitarity is one of the most important criterion 
for any theory to be physically consistent.We start with a brief review of the gauge 


theory of gravity whose gauge symmetries are volume preserving diffeomorphisms. 
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5.1.2 Gauge Theory of the Group of Volume Preserving Dif- 


feomorphisms 


Given our difficulties with gravity, it is important to understand what are sharply 
formulated questions. One might naively expect that a theory of quantum gravity 


should calculate correlation functions — as in QFT - such as 


(O(2)O(Y) ++), Guv(®))  (O(@)O(Y)guv(Z)-++), ete. (5.1) 


where ¢ is some matter field, and g,, is the metric. However, such correlators cannot 
be gauge invariant, and hence physical. The reason is that the gauge symmetries 
of gravity, at least at the semiclassical level, are diffeomorphisms, whose 


infinitesimal forms act on local fields by 


db(x) = E°0,0(2), O09 u(z) = VE) (z), ete. (5.2) 


In field theory, local observables are given in terms of the fundamental fields; but in 
gravity, by this reasoning, no local observables can be gauge invariant. To explore 
an alternative route to a viable field theory of gravitation let us go back to the very 
fundament - namely the two experimental observations that (1) the inertial and grav- 
itational masses of a physical body are numerically equal, mrp = mg, and that (2) the 
inertial energy momentum of a closed physical system is conserved, pf = conserved. 


Taking both together we then can write in the rest frame of the body 


alt 
# = (my,0) 2 (me,0) = ph (5.3) 


where we have tentatively introduced the gravitational energy-momentum pé which 
we keep as an entity a priorily different from p/. Note that p% is conserved due to 
Eqn.(5.3). 
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In GR m; = mq is interpreted as an essential identity which leads to the usual geo- 


metrical description of gravitation. 


We follow a different route and investigate the consequences of viewing m; and mg 
or p and p% as different by their very natures - the prevailing view before Einstein 
which comes at the price of accepting the observed numerical equality my = mq as 


accidential. 


Both pf and p% are four-vectors which are conserved, but in our approach through two 
different mechanisms. Obviously the conservation of pf is related to translation invari- 
ance in space-time. Let us use Noether’s theorem to separately derive the conservation 
of a new four-vector in a field theoretical framework relating it to an inner continous 
symmetry of the theory specified below as invariance under volume-preserving diffeo- 
morphisms. That four-vector is then interpreted as the gravitational energy-momentum 
Pa: 


To mathematically implement the idea we start with infinitesimal diffeomorphisms 


XPS KMS ROE (X) aH 0, 1:9°3 (5.4) 


acting on R‘ with points labelled by X°.X'°(X) denotes an infinitesimal invertible 
and differentiable coordinate transformation of R* and DIF FR‘ acts as a group on 


this space under composition. 


To represent this group on fields we have to add the necessary degrees of freedom to 
the fields in question. So let us take a general matter field @ which has to be defined 
for our purposes on the product of the usual fourdimensional Minkowski space-time 


M? and the four-dimensional inner space R* introduced above 


o(x) > o(@, X) 


with Lagrangian 


Liyg= / d'XA~“*Lu (¢, Ono) (5.5) 


Here, we have introduced a length scale A which is necessary to define an appropriate 


summation over the field components continuosly labelled by X. This summation mea- 
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sure has to be dimensionless as is the summation over the field components discretely 
labelled in the usual case of compact gauge groups - hence the necessity to introduce 
a length scale A so that [ d*XA~ is dimensionless. 


To next represent the group DIF FR‘ on the matter field we take the passive view of 
diffeomorphisms 
pia eee Py SaaS XO (5.6) 
p(x, X) — b'(a, X) = o(x, X) _ E“(X) : VaG(x, X) 


transforming fields, but not coordinates. Defining the variation 6.-. = ..’! —.. of an 


expression under an infinitesimal diffeomorphism or gauge transformation we find 


tes / dEXA'V, - (E*(X) Lar (6, 0,0) = 0 (5.7) 


provided that the infinitesimal gauge parameter €°(X) is divergence-free 
Vob*(X) =0 


This condition reduces the gauge group to volume-preserving diffeomorphisms DI F FR* 
which we will work with in the sequel. It is easily shown to be a group and restricts 
the €°(X) to be divergence-free. 


For the Lagrangian Eqn.(5.5) Noether’s theorem then yields for any solution of the 


field equation four conserved currents 


v = 4 —4 OLM 
Jal) = i aie TC kd (5.8) 


Od (ey =O, a= 0150)3 





and four time-independent charges which is the looked-for conserved inner four-vector 
generating global inner coordinate transformations in field space. After gauging the 
current Eqn.(5.8) will serve as the source of the corresponding gauge field. Note that the 
above construction is transferable to any other given matter field. For clarity a compar- 
ison with the usual Yang-Mills situation might help. In our approach the theory being 


invariant under volume-preserving diffeomorphisms comes along with gravitational mo- 
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mentum conservation - accounted for by field theory means and Noether’s law in anal- 
ogy to how the conservation of e.g. colour in QCD is accounted for through an SU(3) 
global symmetry. In our case, however, the fields live on an infinitedimensional space 
labelled by a continous index X and not as in the QCD case on a three-dimensional 
vector space labelled by a discrete index.The global inner symmetry yields through 
Noether’s theorem the conservation of gravitational momentum, gauging it leads to 
the gauge fields to which we turn next To obtain these we next allow for x-dependent 
infinitesimal volumepreserving gauge parameters €°(X) > E°% (2, X), VaE* (x, X) = 0. 
Note that Eqns.(5.6) still define the representation of the gauge group in field space. 


To ensure the local gauge invariance of the Lagrangian Eqn.(5.5) we have to replace 


ordinary derivatives by covariant ones 0, D, 
DG) =O. FAXES (5.9) 
thereby introducing gauge fields A,,“(x, X) which are divergence-free 


VoA,"(2,X) =0 (5.10) 


consistent with VgE°(x,X) = 0. The transformation law for the gauge field is easily 
derived 
be Ay® = OyE* + Ay? » VgE* — €° - Vp Ay° (5.11) 


and respects V,6-A,,° = 0. By construction we find the Lagrangian Eqn.(5.5) 
FeLae(d,Dyd) == f dXA*V 9 (E%(e,X) Lar (6,Dy6)) =0 (5.12) 


to be invariant under local gauge transformations. The field strength components are 


defined as usual by 


DieX DG NS PV (5.13) 
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They are expressed by the gauge fields as 


Pi yk) = OA, GA) = OA BX) 


(5.14) 
+ AyP(a, X)-V~pAv* (2, X) — Av? (2, X) - Ve Au*(a, X) 
and transform covariantly under a local gauge transformation 
be Pi = Fah = Vee" — 6" = Vek. (5.15) 


Finally, a gauge-invariant, minimal and renormalizable Lagrangian for the gauge field 


is given by 


1 a Vv 
= Taal (x, X) . fe Pe, a) 


Le= / EXK*Le(A 07 AV eA) 


LAA OAL AS 
(5.16) 


5.1.3 Covariant approach and the S-Matrix 


A starting point for discussing the gravitational S-matrix is the observation that, to a 
very good approximation, Minkowski space looks like a good solution of the theory of 
quantum gravity. Indeed, the curvature radius of the observed Universe is of order 10° 
in Planck units - flat space is an excellent approximation for all but the longest-scale 
questions, and we therefore assume that it makes sense to consider the exactly flat 


limit. 


Next, we also observe that there are excitations about this which we describe as ” par- 
ticle” states, such as electrons, photons, etc., whatever their more fundamental de- 
scription may be. Their asymptotic states are described by their momenta p;, and 
other quantum numbers. Moreover, we can consider asymptotic multi-particle states, 
consisting of widely separated particles. And, in a quantum theory, we can ask for the 
amplitudes to transition from a given such ”in” state, in the far past, to another ” out” 


state, in the far future. This collection of amplitudes, which we might denote 


S (pir, Di) = out (pir | Di)in ’ (ook?) 
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is the S-matrix. We should note there are two technicalities that we are glossing 
over. First, for gravity in D = 4 spacetime dimensions, these amplitudes aren’t well- 
defined: instead we must perform a sum over soft gravitons, like with soft photons in 
QED. However, as we will see, this technicality can be avoided in higher dimensions. 
Secondly, a more precise definition of the S-matrix involves a map between normalizable 
Hilbert-space states; momentum-space states give a singlular limit of these, though one 


whose handling is well-understood in asymptotically flat spacetimes. 


Such an S-matrix, which is tightly constrained by properties such as unitarity and 
analyticity, can be a very powerful way to summarize our ignorance of a theory. The S- 
matrix approach is so powerful that it lead to the discovery of string theory: Veneziano 
guessed amplitudes satisfying additional ” duality” properties, and they were then found 
to describe excitations of one-dimensional objects. We might anticipate that such 
study in the context of gravity, supplemented by additional physical input, could bear 


important fruit. 


There are multiple proposed approaches to quantizing gravity, and the S-matrix also 
provides a very important test for these: any theory of quantum gravity should give us a 
means to calculate S, at least approximately, or must provide us with some alternative 


physical quantities. 


This test becomes particularly incisive in the ultraplanckian regime, at center of mass 
energies EF >> Mp, where Mp denotes the D-dimensional Planck mass. That any the- 
ory of quantum gravity must describe this regime follows from very general principles, 
and avoiding this regime appears to require radical assumptions. The first principle is 
that of Lorentz invariance: the Minkowski-space solution is invariant under boosts, and 
we can perform an arbitrarily large boost on a single-particle excitation of it to get a 
state with an arbitrarily large energy. The second principle is a very weak notion of lo- 
cality: one can consider independently boosted particles at very large separations, and 
so prepare a state with a large center of mass energy — and the then allow the particles 
to collide. Plausibly during the big bang particles reached such energies, and it is inter- 
esting to ask whether such collisions are likely to take place anywhere in nature today. 
We observe that cosmic ray accelerators (argued to be active galactic nuclei) produce 
particles up to 1012 GeV, and now and then they could collide; that is only short by 
O(107 ). Even more interestingly, in scenarios with large and/or strongly warped extra 


dimensions, the fundamental Planck mass MD could be as low as the TeV range, and 
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in the most optimistic scenarios black holes could give important signatures at LHC 
While we won’t talk about this phenomenology, there are multiple reviews, including. 
Whether or not this physics is experimentally accessible, the gedanken experiment and 
theoretical problem remain: we need a calculational framework that makes predictions 
in the ultraplanckian regime. As I will describe, here an apparently critical and funda- 
mental problem is unitarity. However, as viewed through the lens of earlier attempts 
to quantize gravity, a more basic concern is that of nonrenormalizability. We will begin 
to understand this problem, and then explore the apparently more fundamental issue 


of unitarity, starting with the perturbative approach to gravity 


5.1.3.1 Perturbative quantization and the Born amplitude 


A starting point for perturbative quantization is the action, 


1 


where FR is the Ricci scalar, Ly, the matter lagrangian, and Newton’s constant Gp is 


related to Mp by the Particle Data Book convention: 
Mp? = (2n)?-*/ (8nGp). 


For simplicity we will typically consider for matter a minimally-coupled scalar of mass 
™m, 


oan 2 1 252 
Lm = —5 (V8)? — 5m? (5.19) 


We investigate perturbations about Minkowski space, 


Suv = Nw + Vf 32nG phy 


and expand in powers of h. This gives an action of the form 


1 a 
S=Smot / d?x {ph (Che + /8Gp [AY Tie? + £L3(h)] +--- \ (5.20) 


Here, quantities evaluated with g = 7 are labeled with zeros, £ is a second-order 
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differential operator (a tensor analog of the laplacian), indices are raised and lowered 
with 7,7 is the stress tensor, and £3 is the third-order term in h. The ellipses contain 
higher powers of @ and h. The propagators are given by inverting the kinetic operators 
for ¢ and h; the latter requires gauge fixing. There are two approaches to this, both 
using the quantity 

hig = lia = Shek (5.21) 


The first is to apply the gauge fixing condition 


Py = 0 (5.22) 


Alternately, one may add a gauge fixing term 


1 > \2 
eae ees (hw) 


to the lagrangian; particularly nice simplifications occur for a = —1. In either case, 
one straightforwardly finds the propagator, and using that can immediately calculate a 
first amplitude, that for tree-level scattering of two matter particles. This is the Born 
amplitude. It is most conveniently given in terms of the Mandelstam invariants; let 
pi, p2 denote the momenta of the incoming particles, and —p3, —p, those for outgoing, 
and define 


s=—(pi + po)? = E’,t =—(ps +p)" = -¢,u= — (pati). 


Here EF is the total center-of-mass (CM) energy, g is the ”momentum transfer,” and 


one finds s +t+u=4m?. If we then define the reduced transition matrix element by 


S=1+iT(s,t)(2m)?6? (~~ P) (5.23) 
and take the limit of large s = E”, we find the high-energy approximation of the Born 


amplitude To : 
To(s,t) ~ —87G ps? /t. 
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5.1.4 Canonical approach and the S-Matrix 


5.1.4.1 Asymptotic States: Quantization of Gauge Field 


In this section we quantize the free volume-preserving diffeomorphism group gauge 
field. Solving the linearized field equations we express all: the fields, the inertial 
energy-momentum and the inner momentum constructed in section 5.1.2 in terms of 
creation and annihilation operators and appropriate polarization vectors. The field 
quanta carry both inertial energy-momentum and inner momentum. We then construct 
the gauge field Fock space and the propagator. Finally we establish the ” gravitational 
limit” to connect the field quanta with observable gauge particles. Let us start with 
the Lagrangian density Eqn.(5.16) for the gauge field defined on M* x M?* 


Lo (Ay®, O° Ay*; Ve Ay") F(x, X) + FM a2, X) 


a 
~ AA? 


A VY a a Qa 
+ yeu AMalz, X) +d Ay (x, X) — sygAMalt, X) > Ay (x, X) 


(5.24) 


adding both a gauge-fixing and a mass term to be able to canonically quantize and to 
deal with eventual infrared problems. The Euler-Lagrange equation for the gauge field 
is 


(8 — p?) AK(x, X) + (1— 20" (0,4Pa(2,X)) = Ja(a,X), (6.25) 
where both gauge field self-interactions and interactions with matter currents have been 


collected in the current J,“ which by construction is divergence-free in space-time as 


well as in inner space 
ONTO (x, X) =0, Nadir Gey) = 0. (5.26) 


Taking the space-time divergence of Eqn.(5.25) 


we find that with or without interactions present the divergence of the gauge field 0,.A°, 


is a free scalar field with mass we This in turn allows us to define the transversal part 
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of the gauge field as usual 


A a a 
Al*(x,X) = A,°(x,X)+ ao! (0°A,*(x, X)) ONAL (x, X) =0 


To keep calculations simple we set the gauge parameter 


oe 


(5.27) 


in the rest of the section. With this choice of \ and with or without interactions present 


the transversal part of the gauge field fulfils 


(—0? — 11”) Aa) SU ee): 


(5.28) 


To canonically quantize we first calculate the canonically conjugate field momenta 


OLG 
0 (0° A, F(x, X)) 
= AW? Fo” g(a, X) + A~?99”0,A’ a(x, X) 





IT ale, X) = 


which live in the gauge algebra like the gauge fields 


VAG) 0, “VPI (e,X).=0 
Defining the delta function transversal in inner space 


T d K | iK(X—-Y KK, 
\ 1. 





we next impose generalized canonical equal-time commutation relations 


12F 


(5.29) 


(5.30) 


(5.31) 


(5.32) 
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consistent with the conditions Eqns.(5.30) on the gauge fields and their conjugate 


momenta. 


It is straightforward to re-express the canonical commutation relations Eqns.(5.32) in 


terms of the gauge fields and their first time derivative 
[An (é, x; eS) A” s(t, Y; Y)| =0 
lA,"(, x; x), Aga(t, Y; Y)| 


Ay*(t,2)X), g(t, ¥)] = imu’ A?P5%9(X — Y)5*(a - y) 


0 (5.33) 


m1 


indicating that the sixteen A,“ are quantized like a set of scalar fields. Some, however, 
will have negative norm in Hilbert space, a fact we will have to deal with below to 


arrive at a sensible theory. 


Next we solve” the field equation Eqn. (5.25) and the constraint Eqn.(5.26) for A = 1 


by Fourier transforming 


F 7 Bk Gis woe 
An@2X0= | mas | Gaye LL 


= Pal (5.34) 





~ 2 
+ 
a 
Q 
a 


ep(k, E(K,T) {a+(k, 7; K,P)e 


and putting k” on the mass shell 


ke= pw? or we =k? +p. (5.35) 


Above y = 0,1,2,3 and [ = 1,2,3 denote the polarization degrees of freedom, 


k 
Ep(k, 0) = . and. .€,(kyy),. “y= 12,3 (5.36) 
and 
RD Pai9 3 (5.37) 


the corresponding polarization vectors in space-time and inner space. Note that [ = 0 


does not appear due to the transversality V,A,°(x, X) = 0 of the gauge field in inner 
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space. €, and €° obey the following relations 
BPs Eplk, 7 # 0) = 0, Ep(k, 7) ver (k, 7) = iq 


3 
kek? 
> e(k, Ne" (k, 1) = —0? + B 





Ky-E%(K,T)=0, €%(K,T)-€.(K,T’) =7™ 





3 
K,K 
S > EK, T)Es(K,T) = —tae + ae for K? >0 
T=1 


To ensure positivity of the gauge field energy the support of the ais restricted to time- 


and light-like vectors in inner space 
supp (a4(k,7;.K,T)) =R* x (V*(K)UV(K)). 


Inversion yields for a+. 


a(n KD) = ian fae fatxa~ 
ciFEHKX 9, (c(h, 9) Ea( KT )A8( ar, X)(t, a3 X)) 


with an analogous expression for ne The commutation relations for the a4, a are 


then easily obtained 


[as (hy KD) ae (hey AP) 0 
jal ( k,7:K,T), at (4,7 HT] =) 





an(b,y3K,T),a) (bys HP)) = Quon? 
A-?(2r)464(K — H)(27)?6°(k — h) 


and show that states have positive norm only for 7, 7/ = 1,2,3. The canonical energy- 
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momentum tensor is calculated as usual 
Bra = paxx (EG OTA oiX | 


1 Vv a Oo 
— ZI” Fys(, X) - FP a(t, X) 
+0,A° a(x, X) 0” AM" (x, X) 
1 
acne 7 As (2, 


x) 
42 
nh Af (x, X) were) )} 
Next we re-express both the momentum three vector 


p= / dx / d*XA °C Al (a, X)-O'A,*(a, X) 


Next we re-express both the momentum three vector 
p= [ate faxrorar(e,X) -O'A,*(x, X) 
and the energy 
= 5 fas fatxas {0° AP (2, X) -0°A,*(2, X) 
+0' A’, (x, X)-0'A,*(a, X) + pw? APa(a, X)- A,*(x, X)} 


i" obtaining the covariant inertial energy-momentum vector 


ol Bk a, i ee 
i ie De 


y=0 Tr 


in terms of the a,,a 








mone {as (k,73 K,T)al(k, 73K, r)+ he. \ 


Note that the energy p° > 0 is positive if the summation is extended over y = 1, 2,3 
only which is the case for purely transversal states. 
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In an analogous way we next determine the inner momentum tensor 
J3(a) = ii dXA° {F#? (x, X) «Vg Ap*(a, X) 
+0,AP a(x, X) - VgA"*(x, X)} 
and the conserved inner momentum vector 


P; = / Bar i. d*X A889 AP, (2, X) - VgA,%(x, X) 
1 dbk i 
=; | oan: i; (nyt Kad d 


nm fas (ky73 KT)ah. (ky KT) + hee. }. 





All of the above formula are correct again up to the introduction of normal ordering. 


The calculation of the commutators of the inertial energy-momentum and the inner 


momentum with a yields 


[Pp al, (k, Ys kK, r)| = qt kal, (k, oF @ r) 


Pay al (k,7;K, r)| Sag Kal, (bk 1) 
allowing for an easy cross-check of the conservation of both types of momenta 
[H,p,]=90, [H,P.] =0 


We establish the usual particle interpretation starting with a vacuum state |0) with 


(0 | 0) = 1 which is annihilated by the destruction operator a+ 
ax(k,y;.K,T)|0) = 0 
and out of which the creation operator al: generates one-particle states 
a'(k, 7; K,T)|0) = one-particle state 


with definite energy-momentum k, inner momentum K and polarizations y, TP and with 
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normalization 
(Ola, (Ra's HP )al. (9 K,1)| 0) = Que ™ 
A? (2r)*64(K — H)(2r)°5°(k —h). 


Acting on a simultaneous eigenstate |h, H) of p, and P, with eigenvalues h,, and Ha 


respectively we find 


pal (ky; K,T)|h, H) =" (hy, + ky) oh (ky; K,T)|h, A) 
Pua+(k, 7; K,T)|h, H) = 777" (hy — ky) 4 (Kk, 7; KT)|h, A) 
Pa, (k,y;K,T)|h, H) =n" (Ha + Ka) a} (k, 5 K,T)|h, H) 
Pas (k, 7; K,T)|h, H) =" (Ha — Ka) 04(k, 7; K,T)|h, A) 














which allows us to construct the Fock space of stationary states by applying multiple 
creation operators on the vacuum state, however, at the price of having to deal with an 
indefinite metric at this point as is usual in the Stueckelberg-Gupta-Bleuler approach 
for quantizing gauge fields. Finally we calculate the time-ordered product of two field 


operators to obtain the Feynman propagator for free gauge fields 


(0|T (Ay*(x, X) A," (2, X))| 0) = —iA?T°? (X —Y). 


i d*k e tk(a—-y) Tw 
(any4 k2 — 2 Lie 





which is local in inner space. Before turning to discuss the ” gravitational limit” we 


define the Fock sub-space of positive norm states requiring their annihilation by areas a 
Fp = {|a) |a° AD )*(x, X)| a) =0}, 


where Ay? is given by 
3 


eas Pk ld ee 
A"@3)= | Gama, | Gay! 


y=0 T=1 





Bpheye ths ay kad jest" ee 
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and its divergence by 
3 


d°k dk 
p A(—)e = 4 
A= | aaa | at aH 


£%(K, TD) (k 0; K, Det 48 





containing only the destruction operators a,for y = 0. Due to Eqn.(101) Fp contains 


positive norm states only 
al (k, 7; K,T)|0) = positive norm states for y,T = 1,2,3 


as required for a consistent probabilistic interpretation of the theory. Finally we intro- 


duce the ” gravitational limit” by mapping K onto k 
al T) = lim al (k, 9; K,T 
a'(k, 7,0) = lim ay(k,7; 4,1). 
The normalization of these states involves again a regularization of the inner volume 


[a(k, 7,0), a4 (B,9/sP)) = Jim Tim Quynh 


= A*(2n)'6"(K — H)(2n)39°(e — hy) 
if / Vy’ 
=Qwogn” aah (2m )°O(k — h) 
=2u, 9" A? (2)°5?(k — h) 
and yields properly normalized destruction and creation operators belonging to a free 


gauge field (but with field quanta still remembering that they carry both inertial and 


inner energy-momentum). 


The physical Fock space Fypys of positive norm states is again obtained requiring their 
annihilation by arASe 


Fowys = {18) |°A,*(2)| 8) = OF, 


where 
3 


Ame) = free DD colle NEE Dab. 7 Te 


3 
y=0 T=1 
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and 
p A(-)a = a ike 
6} Ay (x) lo ae as nye E(k. 1 alk 0, Dye 


The positive norm one-particle states 
a'(k,y,I)|0) = physical states for y, = 1,2,3 


represent then the asymptotic field quanta which correspond to the observable gravitons 


in our approach. 


5.1.4.2 S-Matrix 


In this section we first introduce the operator S mapping the asymptotic in-Fock space 
onto the asymptotic out-Fock space and discuss its properties. We then define the S- 
matrix S' in quantum gravity as the gravitational limit of the operator S mapping the 
physical in-Fock space onto the physical out-Fock space. S is shown to be a unitary 
operator on the asymptotic physical Fock space if S is on the asymptotic Fock space. 


Let us define the operator S by its matrix elements 
Shy Bake: = iy Aig ee cout | bay Bin hey Kn 1) 


or 


Shy Hy ike Ks = (ha, Ay3 3. 65 Rim om, [S|k1, Ky; 2.5 Kay Ky, i ) 


Note that we subsume all quantum numbers other than k and K characterizing an 
asymptotic stationary state in the above notation as they add no new features in our 


context. 


S maps the asymptotic in-Fock space Fi, onto the asymptotic out-Fock space Fout 
ied Seas Welles 1h | Chis Hgsicae ipa sth OU 
and is invertible by construction 
(hits Hache out: |S, 4 aN hitless, at | 
Next we list a few of the key properties of S : 
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A) Invariance of the vacuum yields 
(0 in |S = (0 out | = e*(0 in | 
and its uniqueness allows us to set the phase equal to zero so that 
(0 out | = (0 in | = (0]. 
B) Invariance of the one-particle states yields 


(p, Pin|S|p, P in ) = (p, P out |p, P in ) 
= 49, Pin: |p,P in y=. 


as we have 
|p, P in } =| p, P out ) = |p, P) 


C) The operator adjoint to S is defined by 
SB”) Kh its4 hg Roun \ S| high ees eg Fey Outs) 


and S is by construction unitary 


(hy, Ah; ro oe Tis ded in |SS*| k,, Ky; sa hy ea in ) 
= halls his Bag OUw hig Kae ne ha I Out > 
= On; ,Hyski,Ki 


or 


ssi=sits=1 


if the sets of ortho-normal vectors spanning the asymptotic in- and outFock spaces 
Rigid ha, Kyi) | ki, Ky;...;kn, Kn out ) 


in Eqn. (123) are complete. Note that in the case of gauge fields S is pseudo-unitary 
only and its unitarity on the asymptotic 7n - and out-Fock subspaces of positive norm 
states has to be proven. In the case of the gauge field A,,“(x, X) we have established 
that proof in [10]. 


Next we define the S-matrix S in quantum gravity by the gravitational limit of the 
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matrix elements of the operator S 
Shiki = lim lim Sh; .HyikiKi 


jh; Ki ok; 


= diy lim (hy, gs lig Hy OU Hig Fg ee 5h hy ) 
Hj—-h; Kok; 


= hia t ty oOUr | ieee kins 


The operator S maps the asymptotic physical in-Fock space onto the asymptotic phys- 
ical out-Fock space 
(hy... Nm in |S = (hy... Am out | 


and is invertible 
Cis tig, Cub | SP hy. ep an 


The properties A), B) and C) hold again: A’) Eqn.(118) yields 
(0 in |S = (0 out | = e’7(0 in | 

23 and we can again set the phase equal to zero so that 

0 outs|- == 40 "in’ | 40), 
B’) Eqn.(120) yields 

(p in |S|p in ) = (p out |p in ) 
=(p in |p in )= 1 

as we have again 

| p in ) =| p out ) = |p). 


C’) To prove the unitarity of S on the asymptotic physical Fock space we introduce the 
self-adjoint projection operator P = P? which maps any in- or out-state in the asymp- 
totic positive-norm Fock space on its gravitational-limit state vector in the asymptotic 


physical Fock space Fpnys 


P | Ripka ee hae i, in) = lim le ccee Gp aaeers rm .0n nm). 
iG 
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Note that P is the identity in the asymptotic physical Fock space 
P= Digs 
We can now express S in terms of S and P as 
tes th (PS = cc im |SP. 
Assuming S to be unitary allows us to write 


Lsiys = PIP’ =PS'SP 
= SiPPS = Sts 


This proves the unitarity of the gravitational S-matrix S if S is unitary which is a key 


pre-requisite for a viable physical theory. 


5.1.5 Path Integral Approach to Quantum Gravity 


The basic idea of the Feynman path integral is that the amplitude to go from a state 
with metric g,, and matter fields ¢; at time t, to a state with metric gz and matter 
fields @2 at time tz is given by an integral over all field configurations which take the 


given values at times t, and fo : 


(niet ee = i DlglDId] exp(éT|g, 41) (5.38) 


where D|g] is a measure on the space of all metrics, D[@] is a measure on the space of 
all matter fields, J is the action, and the integral is taken over all field configurations 
with the given initial and final values. (we are using units in which c=h =k = 1. ) 


The gravitational contribution to the action is normally taken to be 


1 


oe i. R(—9)'2d4e (5.39) 


However, the Ricci scalar R contains second derivatives of the metric. In order to 


obtain an action which depends only on first derivatives, as is required by the path- 
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integral method, one has to remove the second derivatives by integrating by parts. 


This produces a surface term which can be written in the form 


1 


K 1/233 
=a} (h)?Ba +C 


where the integral is taken over the boundary of the region for which the action is 
being evaluated, K is the trace of the second fundamental form of the boundary in 
the metric g, h is the induced metric on the boundary, and C’ is a term which depends 
only on the boundary and not on the particular metric g.In order to make sure that 
one registers this surface term correctly one has to join the initial and final spacelike 
surfaces by a timelike tube at some large radius ro. It is convenient to rotate the time 
interval on this timelike tube between the two surfaces into the complex plane so that it 
becomes purely imaginary. This makes the metric on the boundary positive definite so 
that the path integral can be taken over all positive-definite metrics g that induce the 
given metric for the boundary. Suppose that one wants to find the number n(£)dE of 
states of the gravitational and matter fields which have energy between F and E+ dE 
as measured from infinity. This will be given by 


n=— | 2(8)exp(GE)dB (5.40) 


where 


Z(8) = S—(Gns dnl exp(—BH)| gn; bn) (5.41) 


is the partition function for the system consisting of the gravitational and matter fields 
contained in a box of radius 779 at a temperature T = 3~!. This partition function can 
be expressed as a path integral over all matter and gravitational fields that are periodic 
in imaginary time with period {, i.e., 


A 


Z= | Dig)DIelexp(—2) (5.42) 


where [ = —iI is the Euclidean action and the path integral is taken over all positive- 
definite metrics g whose boundary is a two sphere of radius 79 times a circle of cir- 
cumference 3 representing the periodically identified imaginary time axis. One would 


expect that the dominant contribution in the path integral for Z would come from 
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metrics g and matter fields ¢ that are near background fields go,¢ that extremize 
the action, i.e., are solutions of the classical field equations with the given periodicity 
and boundary conditions. Neglecting, for the moment, the question of the radius of 
convergence, one can expand the action in a Taylor series about the background fields 
Ig, ¢| = I (go, go] + L2[9, é]+ higher-order terms where g = gy) + 9,¢ = 60 + ¢, and Ip 
is quadratic in the perturbations g and ¢. If one neglects the higher-order terms, then 


InZ = —f go, da] +10 f Dlg, exp (—Z21¢.4)) (5.48) 


One can regard the first term in the equation above as the contribution of the back- 
ground field to the partition function and the second term as the contribution of thermal 


gravitons and matter quanta the background geometry. 


5.2 Path Integrals and Introductory String Theory 


5.2.1 Sum over worldsheets and The Polyakov path integral 


The Feynman path integral is a natural approach to describe interactions in string 
theory and one way to express a quantum theory. Amplitudes are calculated by sum- 
ming all possible histories and interpolating between the start and end states in route 


integral quantization. Each piece of history is given a certain amount of weight. 


exp (7Sq/h) (5.44) 


with Sq For the given history, the classical action. As seen in figure 5.1(a)for the open 
string and figure 5.1(b) for the closed string, an amplitude is defined in string theory 
by summing over all world-sheets linking specified beginning and end curves. The free 


propagator is produced by adding the equivalent sums for relativistic point particles. 


Strings can interact in a variety of ways, as might be imagined. One example is a 
contact interaction, which occurs when two strings cross. A long-range force mediated 
by a quantum field is another possibility. However, such interactions cannot be added 


to string theory in a fashion that is compatible with the symmetries; we shall see why 
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this is true later. Instead, only interactions that are already implied in the total over 
world-sheets are permitted. Take, for instance, the world-sheets depicted in Figure 5.2. 
Figure 5.2(a) appears to be a quantum correction to figure 5.1(a)’s string amplitude, 
involving intermediate states with two open strings. There are three external variables 
in figure 5.2(b) 


(a) (b) 


Figure 5.1: (a) An open string world-sheet with the topology of a strip. 
The heavier curves are the world-lines of string endpoints. (b) A closed 
string world-sheet with the topology of a cylinder. 


(a) 





Figure 5.2: (a) Quantum correction to open string propagation. (b) Decay 
of one closed string into two. The dashed lines are time-slice. 


closed strings and represents one string decaying into two, it turns out, is the proper 
method for introducing interactions into string theory. We’ll show how these inter- 
actions result in a finite and unitary theory that incorporates gravity. Consider the 
procedure depicted in figure 5.2(b) as observed at different periods. A closed string 


breaks into two, or two strings come together in the reverse procedure. This is the 
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fundamental interaction between closed strings. All particles are produced as various 
states of excitation of the string in closed string theory, and all interactions (gauge, 


gravitational, etc.) are obtained as various states of excitation of the string. 
Polyakov Path Integral 


We now begin to develop the sum over world-sheets, integrating over the fields of the 
Polyakov formalism. We make one change from the Minkowskian world-sheet metric 
Yap is replaced with a Euclidean worldsheet metric gap (o',o7), of signature (+, +). 
The integral runs over all Euclidean metrics and over all embeddings X¥“ (o', 7) of the 


world-sheet in Minkowski spacetime: 


[laxag\exo(-s) 
The Euclidean action is 
S=Sx+2rx 
with 
1 


~ dra’ 
1 Bs 0 1 

y= — | dog!" R+— dsk 
Ar M 21 OM 


The geodesic curvature is denoted by k note that the boundary now is always spacelike. 





Xx ie dag’? og? O,X"O,X, 


(5.45) 


The advantage of the Euclidean path integral is that the integral over metrics is better 
defined. The topologically nontrivial world-sheets we have been describing can have 
nonsingular Euclidean metrics but not Minkowskian ones. We will take the Euclidean 
theory as our starting point, and we will show how to give it a precise definition and 
that it defines a consistent spacetime theory with the advertised properties. However, 
let us give a brief formal argument that it is equivalent to the Minkowski theory with 


which we began. 


We start again with the example of the point particle, where the path integral 


[ianaxtesn | f ar (r'x"X,— nm) 


is oscillatory. The path integral over 7 and X¥ is a product of ordinary integrals, and 
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so we can deform contours just as in ordinary integration. If we take 
n(t) +e n(r), X°(r) + e* X%(r) 


for infinitesimal 6, all terms in the exponent acquire a negative real part and so this 
acts as a convergence factor. Now we can perform a contour rotation in field space, all 


the way to 7 = —ie, X° = —iX”. The integral becomes 


/ [de dX] exp -5 / ie Gs xox rent) 


w=1 


This is the Euclidean point-particle analog of the path integral (3.2.1). We have just 
made a contour rotation, so the Euclidean path integral gives the same amplitude as 


the Minkowski one. 


5.2.2 Path Integrals of Super-symmetric String theory 


Finding a coherent mathematical framework that unifies quantum field theory with 
general relativity is one of the most important concerns in theoretical physics today. 
While the other fundamental interactions, such as electromagnetic, weak, and strong 
interactions, are well characterised by local quantum field theories, gravity has eluded 
any coherent quantum mechanical explanation. At the moment, superstring theory 
is the most promising approach to solving this problem. String theory began as an 
S-matrix theory of hadrons , but it was subsequently revealed that its spectrum always 
contains a massless spin 2 particle, coupling in the same manner as the graviton does 
in general relativity . The bosonic string, which depicts particles with just integer 
spin, is the earliest model. Later, space-time fermions were included into fermionic 
string theory . However, both theories feature a tachyonic particle in their spectra. 
The GSO architecture was used to address this challenge, resulting in a string theory 
with a space-time supersymmetric particle spectrum. This superstring theory quickly 
established itself as a viable contender for a consistent quantum theory that includes 
gravity. String theorists commonly believe that what we now know about superstrings 
is merely a small portion of a much larger structure that has yet to be discovered. It 
is impossible to do any non-perturbative computations or even provide a specification 


of the theory’s ground state in its current form. Superstring theory is still distant 
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from providing predictions that can be tested by experiments because to this poor 
understanding. So far, all hints about what superstring theory is have come from its 
perturbation expansion. String perturbation theory is expressed as an increase in the 
number of loops, much like conventional quantum field theory. Each term is defined 
as a sum over world surfaces with a fixed number of handles, which represent the tra- 
jectories of strings propagating across space-time . A two-dimensional supersymmetric 
quan turn field theory linked in a Weyl invariant way to two-dimensional supergrav- 
ity describes the sum over world surfaces in superstring theory. The interaction of the 
two-dimensional quantum field theory with the surface geometry is finite dimensions 
due to this Weyl invariance. The surface’s non-trivial variations are parametrized by 
a finite collection of coordinates known as superinoduli. Over this supermoduli space, 
string amplitudes are integrals . The representation of Feynman diagrams as sums 
over world lines is the point particle analogue of this description. The string represents 
the supermoduli.Supermoduli are a string generalisation of proper time variables that 
parametrize the length of world lines. A key distinction between strings and point 
particles is that, although point particle interactions have specific locations in space- 
time, string theory interactions are entirely topological; there are no interaction sites 
on a world.surface. Because of this discrepancy, it is predicted that, unlike in quantum 
field theory, the loop corrections for superstring scattering amplitudes will not undergo 
renormalization and will be finite. The geometry of the target space-time is given a 
priori in the path-integral formulation of string theory, and the strings themselves are 
thought of as a disturbance on this classical background. String theory should, in the- 
ory, be able to determine the space-time geometry, thus this is plainly an incomplete 
explanation. Nonetheless, it is an essential topic whether superstring scattering ampli- 
tudes can be described consistently, that is, without infinities, tadpoles, or any other 


form of anomaly, inside a given spacetime. 


A string propagating through space-time sweeps out a two-dimensional surface. For 
bosonic strings the location of this world surface is described by d scalar ficlds x (£1, €2), 
defined on a two-dimensional parameter space. Here we will consider only closed string 
theories, so the parameter space has the form of a cylinder. The superstring has in ad- 
dition to the bosonic coordinates x“ (£1, €2) also intcrnal spin- ning degrees of freedom, 
represented by d fermionic fields 7" (£1, £2). The bosonic and fermionic coordinates 
together constitute a two-dimensional supersymmetric system. By introducing in ad- 


dition a two-dimensional metric g = gnyd&“dEé’ and a gravitino field y,, the sum over 
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fermionic surfaces can be defined as the integral| 6] 


ye = f (agdxaerav e S(gx.et be") 


fermionic 


surfaces 


where S(g,x,2,W) is the action of the (x“,w") -system coupled to two-dimensional 
supergravity [8]. This defiuition of the sum over fermionic surfaces is made invariant 
under reparametrizations and local supersymmetry transformations on the world sheet. 
As a consequence, the functional integral (1.1) has a large redundancy, which has to 


be removed by gauge fixing. We choose the so-called conformal gauge 


Jab = POab Xa = Nake 


In this gauge the action S(g, y,x,w) has the following form: 
1 = 2 
Sof = [ee (5a.0"dat, + UW Oxy + i055) 


Here (z,Z) = ( + t&, € — i€2) are complex coordinate such that g = pdzdz. ’The 
decoupling of the factors p and ¢ from this action is a consequence of the super-Weyl 
invariance of S(g,y,2,~w)[8]. This invariance persists at the qt antum level only in the 
critical dimension d = 10, [6]. This means that, for d = 10, the (a, ~") quantum field 
theory only feels the conformal class of the metric and the gravitino. These conformal 
equivalence classes are called super-Riemann surfaces . We will limit our discussion 
here as the Construction of S-Matrix for any type of String Theory is beyond the scope 
of this thesis. 


144 


CHAPTER 6 


Concluding Remarks and future 
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6.1 Conclusion 


The aim of this work was to make a complete theoretical analysis of the S-Matrix, as 
we are aware that now a days Particle Physicists both experimentalists and Theoreists 
extensively look for new physics coming out of the particle accelerators, where we collide 
a bunch of different particles and observe the zoo of both new and old particles coming 
out of them. It is evident that in most of our collider experiments we have a initial state 
of particles and finally after collision we have a final state of particles, due to the very 
high energy scale and consequently very less time scale it is impossible to tell exactly 
what the particles are doing during the period of collision. As an observer we only get to 
see the final outcoming particles and initial particles that we accelerated.This is where 
the mathematical machinery of S-Matrix comes into play.Computation of S-Matrix 
requires only information about the Initial state of particles prepared in long past and 
final state of particles. We use the asymptotically free approximation and Put several 
condition on our S-Matrix such as Unitarity and Analiticity. These approximations and 
restrictions are enough to Compute the S-Matrix theoretically and study the Scattering 


and collison experiments. 


This is why we developed the theoretical understanding of the S-Matrix in several 
different formulations, fields and contexts. Our discussion started with a very basic 


and brief comparisons of the two major formulations of modern QFT ie the Canonical 
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and Path integral formulation, this was done in order to provide a context as to where 


we will start studying the S-Matrix. 


We then moved towards the study of previous works done by the likes of Dirac, Dyson, 
Weinberg, Schwinger, Tomonaga,etc where we studied the Canonical formulation of the 
S-Matrix.We discussed the S-Matrix as a limit of Time evolution operator and discussed 
its perturbative form also known as the Dyson series. We also studied the properties 
and symmetries of the S-Matrix and finally its Unitarity, Implications and Proof of 


Unitarity towards the end of the chapter. 


We did a complete study of the Path Integral formulation, and discussed the con- 
struction of the Functional form of S-Matrix. AT last we dedicated a whole Chapter 
towards studying various formulations and approaches to the problem of Quantizing 
gravity and constructed the S-matrix of Quantum Gravity. Again this is important 
because S-Matrix, its properties and restrictions may turn out to be a guiding path 
towards a fully consistent theory of Quantum Gravity. We also did a brief study of 
Path Integrals in the Introductory String theory towards the end. 


6.2 Future scope and S-Matrix bootsrap 


The S-matrix bootstrap is a numerical method that can be used to determine or con- 
strain the scattering amplitudes of particles in quantum field theory using simple prin- 
ciples. Over the past few decades, some physicists have tried to use this technique to 
study different physics theories and phenomena. Researchers at Tel Aviv University, 
Ecole Polytechnique Fédérale de Lausanne (EPFL) and the Perimeter institute and 
ICTP-SAIFR have been developing the S-matrix bootstrap further and trying to ap- 
ply it to different areas of physics. In a recent paper published in Physical Review 
Letters, the team tried to use it to investigate string theory, the renowned physics 
theory representing the fundamental components of the universe as one-dimensional 


(1D) ’strings’, instead of point-like particles. 


”The S-matrix bootstrap is an old idea that was popular in the 60s but lost popularity 
with the emergence of quantum chromodynamics to describe the strong nuclear force,” 
Andrea Guerrieri, one of the researchers who carried out the study, told Phys.org. ” The 


goal of the S-matrix bootstrap is to determine (or constrain) the scattering amplitudes 
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of particles in quantum field theory using just basic principles, such as causality, special 
relativity and the fact that probabilities cannot be larger than 1. The novelty we 
brought to the topic was the proposal of a systematic numerical algorithm to implement 


all these principles.” 


While there are now a few dozens of past papers exploring physical systems using 
different versions of the S-matrix bootstrap, most of these works focus on massive 
particles. Guerrieri and his colleagues, however, were recently also able to use it to 
study scattering of massless particles, such as transverse fluctuations of a color flux 


tube in quantum chromodynamics or in massless pions. 


”Once they realized that the method works for massless particles, they tried it for 
gravitons, which are the mediators of the gravitational force,” Subsequently, they chose 
to impose supersymmetry in 10 spacetime dimensions to reduce the technical hurdles 


of the problem.” 


The recent work carried out by Guerrieri, Penedones and their colleague Pedro Vieira 
applies to any theory of quantum gravity. In their study, however, they specifically 
used the S-matrix bootstrap to examine string theory and tried to determine whether it 
is the only consistent quantum gravity theory. ” Einstein’s theory of General Relativity 
(GR) describes the gravitational interaction at long distances or (equivalently) low 
energy,’ Penedones said. ” For example, it can be used to compute how two gravitons 
of large wavelength (much larger than the Planck length) scatter from each other. 
However, if we decrease the wavelength at some point, we cannot trust the prediction 


of GR. In fact, GR’s prediction will give probabilities greater than 1.” 


To improve the reliability of GR’s predictions, the theory must be corrected at short 
distances/wavelength or, equivalently, at high energy. This correction is typically re- 
ferred to as an ‘ultraviolet (UV) completion’. ” We often say that a theory of quantum 
gravity is a UV completion of GR,”. ”In their recent pape, Guerrieri, Penedones and 
Pedro Vieira they investigated the first correction to the GR predictions in the low 


energy limit. 


String theory is a renowned theory of quantum gravity. Notably, in string theory, 
the alpha parameter described by the researchers can take on a specific set of values. 
In string theory, alpha is related to the tension of the fundamental string (in Planck 


units), . They set out to explore in this paper was to be agnostic about what is the 
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theory of quantum gravity which describes the real world and simply ask what values 


of alpha any such theory could produce.” 


Surprisingly, the researchers found that the alpha values that a theory of gravity de- 
scribing the real world could have, according to S-matrix bootstrap numerical calcu- 
lations, were the exact same values produced by string theory. While these findings 
could have interesting implications, their study had a couple of limitations that may 


affect its validity. 


So far, they have found that the values of alpha allowed by general principles coincide 
with those realized in string theory, This is compatible with a positive answer to the 
question above but it is weak evidence. However, in principle, we can study other 
parameters, let’s call them beta, gamma, etc., representing further sub-leading correc- 
tions to the predictions of GR at low energies. There is an infinity of them that String 
Theory predicts. If we can show that the allowed range of these parameters from the 
S-matrix Bootstrap principles also coincides with String Theory, then the evidence will 


be strong.” 
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